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PREFACE 

The  work  reported  here  was  performed  by  the  Optical  Science 
Laboratory  of  the  Advanced  Concepts  Division,  Environmental  Research 
Institute  of  Michigan  (ERIM).  The  work  was  sponsored  by  the  U.S.  Army 
Research  Office  (ARO)  under  Contract  DAAL03-86-K-Q006.  At  ARO,  the 
technical  monitors  were  Dr.  Bob  Guenther  and  Dr.  David  Skatrud. 

This  final  technical  report  covers  work  performed  from  October  1985 
to  October  1988.  The  principal  investigators  at  ERIM  were  Jack 
Cederquist  and  Anthony  M.  Tai .  Major  contributors  to  this  work  were 
Jack  Cederquist,  Anthony  M.  Tai,  Michael  T.  Elsmann,  James  R.  Fienup, 
and  Mark  A.  Stuff. 
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1.0  PROBLEM  STATEMENT 

Optical  systems  employing  a  single  holographic  optical  element 
(HOE)  have  proven  useful  for  many  applications.  Consequently,  many 
theoretical  results  and  computer  simulation  techniques  are  available 
for  the  study  and  design  of  single  HOE  systems.  One  of  the  results 
from  this  work  has  been  an  understanding  of  the  limitations,  both  in 
applications  and  in  performance,  of  single  HOE  systems.  Recently,  a 
few  researchers  have  studied  systems  using  two  or  more  HOEs.  The 
initial  effect  has  been  to  increase  awareness  of  the  potential  new 
applications  and  performance  Improvements  possible  with  such  systems. 
Of  particular  interest  have  been  applications  involving  the  redirection 
of  light  for  beamforming,  optical  coordinate  transformation,  and 
optical  data  interconnection.  The  HOEs  used  have  generally  been  of  the 
computer-generated  type  to  give  the  needed  flexibility. 

With  the  above  results  as  background,  the  goal  of  this  research 
program  was  to  explore  the  possibilities  of  multiple  HOE  systems  when 
used  for  generalized  beamforming.  Generalized  beamforming  includes 
both  the  redistribution  of  light  Intensities  and  the  more  demanding 
one-to-one  mapping  of  optical  data  coordinate  transformation.  Specific 
goals  were  to  increase  the  theoretical  understanding  of  generalized 
beam  forming  with  multiple  HOEs,  to  develop  practical  methods,  whether 
theoretical  or  algorithmic,  for  determining  the  phase  functions  which 
specify  the  individual  HOEs,  to  investigate  applications  for  such 
systems,  and  to  fabricate  and  experimentally  investigate  the 
performance  of  multiple  HOE  systems  for  generalized  beamforming.  It  is 
understood  that  the  required  HOEs  will  need  to  be  computer-originated, 
that  is,  at  some  stage  of  their  fabrication,  computer-generation 
techniques  will  be  necessary  to  realize  the  complicated  phase 
functions. 
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2.0  SUMMARY  OF  RESULTS 

2.1  HOLOGRAPHIC  POLAR  FORMATTING 

A  two  holographic  optical  element  (HOE)  system  for  polar  formatting 
of  spotlight  mode  synthetic  aperture  radar  (SAR)  data  was  designed, 
fabricated,  and  successfully  tested.  With  the  addition  of  a  spatial 
light  modulator,  a  third  phase-compensating  HOE,  and  a  Fourier 
transform  lens,  the  real-time  polar  formatting  of  SAR  data  and  SAR 
image  formation  was  experimentally  demonstrated. 

SAR  image  formation  is  an  excellent  application  for  coherent, 
analog  optical  processing.  However,  for  spotlight  mode  SAR,  polar 
formatting  of  the  SAR  data  is  required.  One  method  for  real-time 
optical  polar  formatting  would  be  to  use  coordinate  transform  HOEs 
designed  using  the  stationary  phase  approximation.  Unfortunately,  a 
single  HOE  cannot  perform  the  polar  coordinate  transformation  for  the 
SAR  data  processing  application.  However,  it  was  found  that  polar 
formatting  can  be  performed  optically  by  the  successive  use  of  two 
holographic  coordinate  transformations.  The  approach  Investigated  and 
demonstrated  uses  a  recti  1 i nearly  addressed  SLM  to  input  the  SAR  data 
and  a  two  HOE  system  to  optically  polar  format  the  data.  By  using  a 
recti  linearly  addressed  SLM,  this  approach  will  benefit  from  the 
current  intensive  development  of  electronically  addressed  SLMs. 

Briefly,  the  system  operates  as  follows:  The  SAR  data  was  input 
either  by  recording  it  on  film  or  by  using  an  electronically  addressed, 
liquid  crystal  television  spatial  light  modulator.  In  either  case,  the 
data  was  Imaged  onto  the  optical  polar  formatter  consisting  of  two  HOEs 
and  two  Fourier  transform  lenses.  A  third  HOE  Is  necessary  to  correct 
the  output  phase  so  that  a  final  lens  can  perform  a  Fourier 
transformation  to  produce  the  desired  SAR  image.  The  phase  functions 
for  the  HOEs  were  implemented  as  computer-generated  holograms  (CGHs)  by 
encoding  them  as  real,  non-negative  amplitude  transmittance  functions. 
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To  Improve  their  diffraction  efficiency,  the  CGHs  were  copied  onto  a 
volume  phase  material  and  placed  on  a  high  carrier  frequency.  The 
result  is  a  computer-originated  holographic  optical  element  (COHOE). 

The  multiple  HOE  system  was  used  experimentally  to  polar  format  and 
process  SAR  data  from  multiple  point  scatterers  with  both  the  film  and 
the  liquid  crystal  television  input.  The  experimental  results 
demonstrated  high  SNR,  low  aberrations,  and  acceptable  efficiency.  The 
theoretical  space-bandwidth  product  (SBWP)  of  the  system  (given  by  the 
square  root  of  the  CGH  SBWP)  was  achieved  experimentally.  Further 
details  are  given  in  Appendix  A. 

2.2  HOLOGRAPHIC  BEAMFORMER 

A  design  method  was  developed  for  computing  the  phase  functions  of 
an  energy-efficient  system  using  two  holographic  optical  elements 
(HOEs)  for  converting  a  Gaussian  beam  into  a  uniform  beam  with 
rectangular  support  in  the  far-field  of  the  source.  The  method  is 
based  on  a  modification  of  the  Gerchberg-Saxton  algorithm  which 
includes  an  x-y  separability  constraint  on  the  phase  of  one  of  the 
HOEs.  A  beamforming  system  was  fabricated  using  this  method,  and 
experimental  results  were  obtained  which  support  the  design  approach. 

In  applications  such  as  laser  radar  and  active  coherent  imaging,  it 
is  often  desired  to  uniformly  illuminate  a  rectangular  area  which  is  in 
the  far-field  of  the  source.  We  investigated  the  use  of  a  holographic 
beamforming  system  to  provide  the  energy  redistribution  into  a  near- 
field  sine  (x)$1nc  (y)  intensity  profile  with  the  correct  phase.  In 
order  to  provide  both  the  energy  redistribution  and  phase  correction, 
two  separated  HOEs  are  required.  The  beamformer  concept  which  was 
developed  therefore  uses  a  first  HOE  to  redistribute  the  energy  of  the 
input  beam  into  the  necessary  near-field  intensity  profile,  and  then 
uses  a  second  HOE  to  place  the  correct  phase  onto  this  beam.  By 
relying  solely  on  energy  redistribution  and  phase  correction,  the 
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system  Is  limited  In  efficiency  only  by  the  maximum  achievable 
diffraction  efficiency  of  the  HOEs,  which  can  be  very  high. 

After  determining  a  near-field  amplitude  distribution  (a  weighted 
sine)  which  yields  the  desired  far-field  characteristics,  the  design  of 
the  beamforming  system  consists  of  computing  the  required  phases  for 
the  two  HOEs  to  produce  the  desired  near-field  beam.  The  first  design 
method  considered  was  adapted  from  that  of  Bryngdahl  for  producing 
geometric  transformations  with  a  single  HOE.  This  stationary  phase 
design  for  the  first  HOE  was  tested  by  computer  simulation  and  found  to 
be  inadequate  for  accurately  producing  the  zeros  of  the  distribution 
which,  in  turn,  resulted  in  unacceptable  ringing  and  high  sidelobes  in 
the  far-field. 

Because  of  the  Fourier  transform  relationship  between  the  two 
planes,  the  design  problem  is  equivalent  to  that  of  computing  an  object 
phase  for  a  known  object  modulus  (truncated  Gaussian)  and  a  known 
Fourier  modulus  (weighted  sine).  This  synthesis  problem  Is  equivalent 
to  that  of  phase  retrieval  from  a  pair  of  Intensity  measurements,  in 
which  the  complex  fields  In  two  domains  related  by  a  Fourier  transform 
are  reconstructed  from  Incomplete  information  in  each  domain.  A  2-D 
Iterative  synthesis  of  the  first  HOE  phase  using  the  Gerchberg-Saxton 
algorithm  was  therefore  attempted.  The  algorithm  was  modified  to 
enforce  a  separability  constraint  on  the  object  phase.  The  stationary 
phase  design  was  used  as  an  initial  phase  estimate  to  start  the 
Gerchberg-Saxton  algorithm,  and  a  specified  number  of  iterations  were 
performed.  A  1-0  slice  of  the  current  phase  estimate  was  then 
extracted,  and  an  outer  product  performed  to  produce  a  2-D  separable 
phase  estimate.  This  was  used  to  start  another  set  of  iterations,  and 
the  loop  repeated  until  an  acceptable  solution  was  found  by  examining 
the  simulated  far-field  Intensity.  After  each  separability  constraint 
Is  Imposed,  the  phase  estimate  Improves,  so  that  overall  convergence  to 
a  good  separable  phase  solution  is  obtained. 
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To  verify  the  results  of  the  computer  simulation,  the  first  HOE  was 
produced  as  a  COHOE  using  a  C6H  having  the  phase  synthesized  by  the 
modified  Gerchberg-Saxton  algorithm.  The  second  HOE  was  produced  by 
interfering  the  output  of  the  first  element  with  a  beam  carrying  the 
desired  near-field  phase  of  the  sine  distribution.  This  desired  phase 
(alternating  r  phase  shifts)  was  encoded  and  recorded  as  a  CGH.  The 
experimental  results  agreed  with  the  simulation  results  sufficiently 
well  to  support  the  design  method.  Further  details  are  given  in 
Appendix  B. 

2.3  REALIZABLE  COORDINATE  TRANSFORMATIONS 

Optical  processors  often  operate  on  data  that  are  presented  in  two 
spatial  dimensions.  In  some  applications,  the  data  need  to  be 
spatially  rearranged  at  some  point  In  the  processing.  These  data 
rearrangements  may  be  thought  of  as  one-to-one  coordinate 
transformations.  Conventional  optics  cannot  easily  perform  complicated 
transformations,  so  holographic  optical  elements  (HOEs)  are  often  used 
for  coordinate  transformations.  Unfortunately,  many  coordinate 
transformations  of  Interest  In  optical  processing  are  not  realizable 
with  a  single  HOE.  An  example  is  polar  formatting  for  image  formation 
from  spotlight  mode  SAR  data.  However,  two  HOEs  in  series  can  be  used 
to  perform  this  coordinate  transformation.  In  this  case,  the  desired 
transformation  Is  decomposed  Into  two  serial  transformations,  each  of 
which  is  realizable  with  a  single  HOE. 

Given  this  example,  two  questions  naturally  arise:  First,  can  all 
one-to-one  coordinate  transformations  satisfying  some  modest 
assumptions  be  achieved  by  the  use  of  multiple  HOEs  In  series?  Second, 
If  so,  what  Is  the  minimum  number  of  HOEs  required?  It  was  shown  that 
all  coordinate  transformations  which  are  sufficiently  differentiable 
can  be  decomposed  Into  two  transformations  in  series  and  therefore 
performed  using  two  holographic  elements.  Although  the  method 
developed  for  designing  the  two  required  holographic  elements  is  not 
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explicit  rjr  the  general  case,  a  useful  necessary  condition  for  the 
decomposition  was  found.  A  simple  explicit  solution  based  on 
separability  of  the  necessary  equations  for  the  first  of  the  two  serial 
transformations  was  developed  and  is  useful  in  some  cases.  Further 
details  are  given  In  Appendix  C. 

2.4  CGH  FABRICATION  FACILITIES 

Because  of  the  Importance  of  computer-generated  holograms  (CGHs)  to 
this  research,  a  small  part  of  the  total  effort  was  devoted  to  a 
consideration  of  existing  CGH  fabrication  facilities  and  techniques. 
Although  many  facilities  are  currently  fabricating  CGHs  of  various 
types,  all  of  these  facilities  suffer  from  one  or  more  of  the  following 
deficiencies:  (1)  performance  limitations  (e.g.,  resolution,  space- 
bandwidth  product,  or  recording  speed),  (2)  cost,  (3)  CGH  encoding  or 
materials  restrictions  (e.g.,  restriction  to  binary  amplitude  or  phase, 
or  restriction  to  recording  fringes),  (4)  data  format  inefficiencies, 
and  (5)  lack  of  technical  support  In  CGH  design  and  optics  in  general. 
These  deficiencies  raise  the  Issue  of  the  need  for  a  national  facility 
for  the  fabrication  of  CGHs.  Further  details  are  given  in  Appendix  D. 

The  SPIE  Symposium  In  Los  Angeles  In  January  1988  included 
conferences  entitled  Holographic  Optics:  Design  and  Applications  and 
Computer-Generated  Holography  II.  To  benefit  those  attending  these 
conferences,  an  evening  panel  discussion  entitled  "The  Need  for  a 
National  CGH  Facility"  was  organized  and  held.  The  consensus  of  the 
panel  members  was  that  software  difficulties  (e.g.,  data  format 
Inefficiency  and  Incompatibility)  are  the  most  Important  limitations  on 
the  development  and  use  of  CGHs  and  that,  for  the  near  term,  electron 
beam  writers  are  the  preferred  hardware  for  the  fabrication  of  high 
performance  CGHs. 
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HOLOGRAPHIC  POLAR  FORMATTING  AND  REAL-TIME  OPTICAL  PROCESSING 
OF  SYNTHETIC  APERTURE  RADAR  DATA 

(submitted  to  Applied  Optics,  1989) 

Abstract 

A  two  holographic  optical  element  (HOE)  system  for  polar  formatting 
of  spotlight  mode  synthetic  aperture  radar  (SAR)  data  was  designed, 
fabricated,  and  successfully  tested.  With  the  addition  of  a  spatial 
light  modulator,  a  third  phase-compensating  HOE,  and  a  Fourier 
transform  lens,  the  real-time  polar  formatting  of  SAR  data  and  SAR 
image  formation  was  experimentally  demonstrated. 

I.  Introduction 

Synthetic  aperture  radar  (SAR)  image  formation  is  an  excellent 

application  for  coherent,  analog  optical  processing.  However,  for 

spotlight  mode  SAR,  polar  formatting  of  the  SAR  data  Is  required  before 

the  image  can  be  formed  via  2-D  Fourier  transformation.  Previous  work 

has  Investigated  the  use  of  mechanical  rotation  and,  in  one  case, 

spatial  light  modulators  with  electron  beam  writing  to  perform  polar 
12  3 

formatting.  '  '  We  view  mechanical  motion  as  undesirable  and  wish  to 
take  advantage  of  the  current  intensive  development  of  SLMs  with 
rectillnearly  placed  pixels  and  electronic  addressing.  We  therefore 
seek  optical  methods  for  polar  formatting.  Single  holographic  optical 
elements  (HOEs)  to  perform  coordinate  transformations  were  first 
designed  using  the  stationary  phase  approximation  by  Bryngdahl4,5  in 
1974.  Unfortunately,  a  single  HOE  cannot  perform  the  polar  coordinate 
transformation  for  the  SAR  data  processing  application.  However,  the 
use  of  two  coordinate  transform  HOEs  In  series  has  been  investigated6,7 
to  overcome  the  limitations  of  single  HOE  systems. 


In  this  paper,  the  design,  fabrication,  and  testing  of  a  three  HOE 
system  for  the  real-time  polar  formatting  of,  and  image  formation  from, 
spotlight  SAR  data  is  described.  In  Sec.  II,  background  information  on 
SAR  data  processing  is  given.  The  design  of  the  required  HOEs  for 
polar  coordinate  transformation  Is  described  in  Sec.  III.  In  Sec.  IV, 
the  fabrication  of  the  HOEs  and  the  construction  of  the  optical  system 
is  discussed.  Experimental  results  for  both  polar  formatting  and  real¬ 
time  SAR  image  formation  are  given  in  Sec.  V.  Conclusions  and 
suggestions  for  further  research  are  presented  in  Sec.  VI. 

II.  Synthetic  Aperture  Radar  Data  Processing 

In  synthetic  aperture  radar,  an  appropriate  waveform  is  transmitted 
and  the  return  signal  scattered  from  a  scene  is  processed  to  obtain 

O 

range  and  Doppler  information  about  the  scene.  Range  resolution  is 
achieved  in  Inverse  proportion  to  the  frequency  bandwidth  of  the 
transmitted  signal.  The  Doppler  Information  corresponds  to  azimuthal, 
or  cross-range,  spatial  information.  Azimuthal  resolution  is  achieved 
in  inverse  proportion  to  the  time  interval  over  which  return  signals 
are  collected  as  the  SAR  flies  by  the  scene  (i.e.,  in  inverse 
proportion  to  the  extent  of  the  synthetic  aperture).  Together,  the 
range  and  Doppler  Information  give  a  2-D  image  of  the  scene.  For  a  SAR 
which  transmits  a  sequence  of  linear  FM  (chirp)  pulses,  the  required 
data  processing  consists  of  four  steps.  The  first  mixes  the  return 
signals  with  a  reference  signal  to  demodulate  the  signal  to  an 
Intermediate  frequency.  The  return  signal  from  each  pulse  thus  gives  a 
1-D  data  set.  The  second  step  forms  a  2-D  data  set  by  placing  the  1-D 
data  from  each  of  many  pulses  side  by  side.  In  the  Important  cases 
where  the  scene  consists  of  a  rotating  object  or  the  SAR  is  used  in  the 
spotlight  mode,  Fig.  1(a)  shows  the  resulting  data  for  a  single  point 
scatterer.  In  these  cases,  the  third  step  polar  formats  the  data  as 
shown  In  Fig.  1(b)  to  compensate  for  the  continuously  changing  viewing 
angle  during  data  collection.  The  fourth  step  performs  a  2-D  Fourier 
transformation  of  the  data  to  produce  an  image  of  the  scene. 
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The  necessary  polar  formatting  operation  is  a  difficult  data 
processing  step.  In  any  sampled  data  processor  (such  as  an  electronic 
digital  computer),  polar  formatting  requires  Interpolation  which  often 
introduces  error  and  is  computationally  Intensive,  about  twice  that  of 
the  FFT  operation.  Since  a  spatial  light  modulator  (SLM)  must  be  used 
to  input  the  data  to  an  optical  processor,  one  possibility  that  has 
been  investigated  is  to  polar  format  the  data  by  writing  it  onto  a 
rotating  SLM  with  an  electron  beam.1  An  optically  addressed  SLM  could 
also  be  used,  but  it  must  have  a  storage  time  as  long  as  the  time  taken 
to  write  all  the  2-D  data  coupled  with  a  short  erasure  time  to  allow 
real-time  input  of  the  next  2-D  data  set.  Other  possibilities  such  as 

rotating  the  output  detector  of  the  data  processor  and  processing  each 

2  3 

l-D  data  set  separately  have  also  been  Investigated.  ' 

The  approach  investigated  and  demonstrated  in  this  paper,  however, 
uses  a  recti  11  nearly  addressed  2-D  SLM  to  input  the  data  in  the  form 
shown  in  Fig.  1(a)  and  a  two  HOE  system  to  optically  polar  format  the 
data  as  shown  in  Fig.  1(b).  It  does  not  suffer  from  the  frame  rate 
limitations  of  optically  or  electron  beam  addressed  SLMs  and  does  not 
involve  mechanical  motion.  By  using  a  rectillnearly  addressed  SLM, 
this  approach  will  benefit  from  the  current  intensive  development  of 
electronically  addressed  SLMs,  such  as  pixel  addressed  liquid  crystal 
light  valves  and  the  S1-PLZT  device.  These  devices  are  compact, 
have  controllable  storage  and  erasure  rates,  and  are  relatively  fast. 
Together  with  optical  polar  formatting  and  optical  processing,  they 
have  great  potential  for  application  to  real-time  processing  of  SAR 
data. 

III.  Design  of  Polar  Coordinate  Transformation  HOEs 

As  an  introduction  to  the  design  method,  consider  the  optical 
system  shown  in  Fig.  2  where  the  input  field  at  Pj  in  the  front  focal 
plane  of  a  lens  L  is  unlphase  and  has  an  amplitude  distribution  a^x, 
y).  A  HOE  with  phase  0(x,  y)  is  placed  in  plane  Pj.  The  complex 
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amplitude  distribution  a2(u,  v)  found  In  the  back  focal  plane  P2  of  the 
lens  Is  given  approximately  by  the  Fourier  transform: 


a2 Cu f v)  =  “  {?  j"  /"^(x.y)  exp[i>(x,y)  -  ^  (xu  +  yv)]  dxdy  (1) 

where  X  is  the  wavelength  and  f  is  the  lens  focal  length.  Following 
4  5 

Bryngdahl ,  '  this  Integral  can  be  approximated  by  the  method  of 
stationary  phase11  to  show  that  a  point  (x,  y)  in  the  input  is  mapped 
to  a  point  (u ,  v)  in  the  output  where 


[Xf 

x' 

u 

2  r 

9x 

y. 

▼ 

V. 

Xf 

-^(x,y) 

L2r 

9y 

(2) 


If  a  desired  coordinate  transformation  is 


g(x,  y) 

h(x,  y) 


(3) 


then  the  required  HOE  phase  ^(x,y)  Is  found  by  solving  the  set  of 
partial  differential  equations: 


UU  =  a(x  y) 

2t  9x  g^x'  y) 

hltt  .  hfx  y) 

2t  9y  ‘ 

1 2 

Equation  (4)  has  a  solution  If  and  only  If 

9q(x.  v)  _  9h(x,  v) 

9y  "  9x 


(4) 


(5) 
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Unfortunately,  for  many  coordinate  transformations  of  interest,  Eq. 
(5)  Is  not  satisfied.  Specifically,  for  polar  formatting,  the  desired 
coordinate  transformation  is 


x' 

u' 

x  cos  y' 

= 

■y. 

.V. 

.x  sin  y. 

and  the  transformation  functions  are 


g(x,  y)  =  x  cos  y 
h(x,  y)  =  x  sin  y  . 


The  partial  derivatives, 


sg(xt-  v) 

ay 

ah(x,  y) 
ax 


-x  sin  y 
sin  y  , 


(6) 


(7) 


(8) 


do  not  satisfy  Eq.  (5)  and,  hence,  the  transformation  is  not  realizable 
with  a  single  coordinate  transform  HOE. 


However,  the  polar  formatting  transformation  can  be  separated  into 
two  cascaded  transformations,  each  of  which  does  satisfy  Eq.  (5).13'14 
The  transformations  are: 


x1 

'u' 

'In  x' 

-► 

= 

yJ 

V. 

y  ■ 

and 


(9) 


u' 

s' 

exp  u  cos  v' 

= 

,v. 

,t. 

.-exp  u  sin  v. 

(10) 
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Therefore,  both  transformations  are  realizable  optically.  Note  that 
the  minus  sign  in  Eq.  (10)  is  necessary  to  make  the  second 
transformation  realizable.  It  corresponds  to  an  inversion  in  one  of 
the  output  coordinates.  Similar  use  of  a  minus  sign  in  Eq.  (6)  does 
not  make  that  transformation  realizable,  however. 


Polar  formatting  can  therefore  be  performed  optically  by  the 
successive  use  of  two  holographic  coordinate  transformations  as  shown 
in  Fig.  3.  HOEs  HI  and  H2  perform  coordinate  transformations  which, 
when  cascaded,  produce  the  polar  coordinate  transformation.  HOEs  HI' 
and  H2 '  correct  the  phase  at  the  output  of  each  transformation  such 
that  the  output  phase  is  uniform  for  a  uniphase  input.  It  should  be 
noted  that  the  phase  correction  for  the  output  of  the  first  coordinate 
transformation  and  the  phase  for  the  second  transformation  can  be 
combined  and  Implemented  as  a  single  HOE,  so  the  complete  system, 
Including  the  final  phase  correction,  requires  only  three  HOEs. 
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The  next  step  Is  to  solve  for  the  phase  functions  of  the  two 
coordinate  transformation  HOEs.  To  do  this,  we  need  to  define 
coordinate  systems  in  each  of  the  three  optical  planes  in  Fig.  3  and  to 
establish  the  relationships  of  these  systems  to  coordinate  systems  for 
the  SAR  data.  As  explained  in  Sec.  II,  after  the  second  step  in  the 
data  processing,  the  Input  SAR  data  lie  in  a  rectangular  region  from  r1 
to  r 2  and  -0^  to  9^  in  a  coordinate  system  in  which  r  and  9  are 
orthogonal  [see  Fig.  4(a)].  This  SAR  data  will  be  loaded  into  an 
electronically  addressed  SLM  with  rectilinear  pixel  placement  such  that 
the  return  signal  data  from  each  pulse  will  reside  in  one  row  of  SLM 
pixels.  In  the  optical  plane  located  at  the  HOE  HI,  this  data  can  be 
placed  In  a  region  defined  to  be  from  -Xj  to  Xj  and  -yj  to  yj  [again, 
see  Fig.  4(a)].  The  relationship  between  these  two  coordinate  systems 
(one  optical  and  one  for  the  input  SAR  data)  is: 


and 


r2  +  rl1 


r 


r2  +  rj 


9 


V 
*1  * 


(13) 


(14) 


The  first  HOE  HI  will  transform  [see  Eq.  (9)]  the  SAR  data  into  a 
region  from  In  r^  to  In  r2  and  -9^  to  in  a  coordinate  system  in 
which  In  r  and  9  are  orthogonal.  In  the  optical  system,  this  data  can 
be  placed  in  a  region  from  -Uj  to  Uj  and  -Vj  to  v1  [see  Fig.  4(b)]. 
The  relationship  between  these  data  and  optical  coordinate  systems  at 
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the  output  plane  of  the  first  coordinate  transformation  (at  the  plane 
of  the  HOE  incorporating  HI1  and  H2)  is: 


2u. 


u  = 


v  = 


In  (r*/^) 


'  In  (r2rj 

In  r - =2-i- 


V 


(15) 


and 


r  =  Jr2r1  exp 


fu  In  (r2/r1)'| 


2u, 


e  = 


v 


(16) 


Combining  Eqs.  (14)  and  (15)  gives  the  required  form  of  the  first 
optical  coordinate  transformation: 


(17a) 


x' 

u' 

-► 

y. 
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where 
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(17b) 


Note  that  Eq.  (17)  is  simply  a  more  general  form  of  Eq.  (9)  which  still 
satisfies  Eq.  (5).  Solving  the  partial  differential  Eq.  (4)  by  using 
the  expressions  In  Eq.  (17)  and  Integration  gives  the  required 
holographic  phase  function  for  the  HOE  HI: 
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(18) 


The  second  transformation  [see  Eq.  (10)]  will  put  the  SAR  data  into 
a  region  from  ^  to  r2  and  -0^  to  0 ^  in  a  coordinate  system  in  which  r 
and  0  are  polar  coordinates.  In  the  optical  system,  the  data  will  be 
placed  in  a  region  from  -s1  to  Sj  on  the  s-axis  and  -tj  to  tj  on  the 
arc  corresponding  to  r  =  (r1  +  r2)/2  [see  Fig.  4(c)].  The  relationship 
between  these  data  and  optical  coordinate  systems  is: 


2s. 


s  = 


t  = 


r2  -  rl 


r  cos  0  - 


'VU'i1 


-2^  sin  0 
(r2  +  rl]  s1n  9 1  * 


(19) 


It  is  convenient  to  require  that 


(r2  +  rl]  sin  5] 
r2  -  rl 


(20) 


so  that  SAR  image  produced  is  scaled  identically  in  both  coordinates. 
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Combining  Eqs.  (16)  and  (19)  gives 
coordinate  transformation: 
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Note  that  Eq.  (21)  is  simply  a  more  general  form  of  Eq.  (10)  which 
still  satisfies  Eq.  (5)  (when  written  In  terms  of  u-v  coordinates) 
provided  that 


u 

v 


1 

1 


In  (r^rj) 


(22) 


Eq.  (22)  expresses  the  requirement  that  the  SAR  data  be  scaled 
identically  In  both  coordinates  in  the  u-v  plane.  In  these 
experiments,  Uj  was  chosen  equal  to  Vj  and  r2/r1  =  5/3,  so  9 1  =  14.6°. 
Solving  the  partial  differential  Eq.  (4)  (In  u-v  coordinates)  by  using 
the  expressions  given  In  Eqs.  (21)  and  (22)  and  integration  gives  the 
required  holographic  phase  function  for  the  HOE  H2: 
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IV.  HOE  Fabrication  and  Optical  Setup 

The  phase  functions  of  Eqs.  (18)  and  (23)  were  implemented  as 
computer-generated  holograms  (CGHs)  by  encoding  them  as  real,  non¬ 
negative  amplitude  transmittance  functions  t(x,  y)  of  the  form: 

t(x,  y)  =  b  +  a  cos  0(x,  y)  +  ax  +  py]  (24) 

where  a  and  b  are  the  amplitude  and  bias  of  the  cosinusoidal  fringe 
pattern.  For  these  experiments,  a  =  0.4  and  b  =  0.5.  The  term  of  the 
form  ax  +  py  was  added  to  each  phase  to  give  a  spatial  carrier 
frequency,  a  and  p  were  chosen  so  that,  in  a  Fourier  plane  with 
respect  to  the  CGH,  the  desired  first  diffraction  order  (the  coordinate 
transform  region)  did  not  overlap  with  the  zero  order.  The  amplitude 
transmittance  functions  were  recorded  on  high  resolution  film  using  a 
rotating  drum  laser  beam  recorder.  A  CGH  size  of  ~(1  cm)  was  used 
with  a  5  /im  pixel  size  to  produce  a  CGH  space-bandwidth  product  of  2048 
by  2048.  The  film  was  mounted  between  1/4  inch  microflat  plates  using 
a  UV-curing  optical  cement  in  order  to  eliminate  any  phase  errors  due 
to  thickness  variations  In  the  film. 

These  CGHs  are  thin  absorption  holograms  and  are  therefore  too  low 
in  diffraction  efficiency  (-1%)  to  be  used  in  a  multiple  HOE  system. 
To  Improve  their  diffraction  efficiency,  the  CGHs  were  copied  onto  a 
volume  phase  material  and  placed  on  a  high  carrier  frequency.  The 
result  Is  a  computer-originated  holographic  optical  element 
(COHOE) .  *5, The  optical  setup  for  producing  a  C0H0E  is  shown  in  Fig. 
5.  A  laser  beam  is  split  Into  two  parts.  One  part  is  collimated  and 
used  to  read  out  the  CGH.  The  desired  diffracted  order  of  the  CGH  is 
passed  by  a  Fourier  plane  spatial  filter  and  imaged  to  the  C0H0E 
recording  plane.  There  It  Is  Interfered  with  a  reference  plane  wave  at 
an  offset  angle  suitable  for  forming  a  volume  hologram.  When 
illuminated  by  the  reference  plane  wave,  the  phase  distribution  encoded 
in  the  CGH  is  reconstructed  with  high  efficiency  by  the  C0H0E. 
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As  shown  In  Fig.  3,  two  coordinate  transform  HOEs  and  a  third 
phase-compensating  HOE  are  used  in  the  polar  formatting  system.  The 
first  HOE  is  a  COHOE  of  the  phase  function  for  the  first  coordinate 
transformation  [Eq.  (18)].  The  second  HOE  corrects  the  output  phase  of 
the  first  coordinate  transformation  and  produces  the  phase  for  the 
second  transformation  [Eq.  (23)].  This  second  HOE  was  fabricated  by 
interfering  the  output  wavefront  of  the  first  HOE  (for  a  plane  wave 
input)  with  the  desired  phase  produced  by  the  second  CGH.  It  is 
therefore  a  COHOE  made  with  a  nonplanar  reference  wavefront.  The  third 
HOE  corrects  the  output  phase  of  the  second  coordinate  transformation. 
It  was  made  by  interfering  the  output  wavefront  from  the  the  second  HOE 
with  a  uniform  reference  beam.  All  three  HOEs  were  made  on  Kodak  649F 
spectroscopic  plates  using  a  HeNe  laser  (X  =  632.8  nm) .  The  HOEs  were 
processed  by  the  silver  halide  (sensitized)  gelatin  process17  and  then 
covered  with  microflats  to  avoid  degradation  of  diffraction  efficiency 
when  exposed  to  humidity.  The  resulting  diffraction  efficiencies  for 
the  three  HOEs  were  37%,  14%,  and  36%,  respectively.  Because  the 
second  and  third  HOEs  are  made  by  Interfering  beams  of  non-uniform 
intensity,  the  recording  must  be  confined  to  the  linear  exposure  region 
to  avoid  Intermodulation  noise.  This,  of  course,  decreases  the  maximum 
diffraction  efficiency  that  can  be  achieved. 

The  experimental  multiple  HOE  system  for  polar  formatting  and  SAR 
image  formation  Is  shown  In  Fig.  6.  The  SAR  data  was  input  either  by 
recording  It  on  film  and  placing  It  in  a  liquid  gate  or  by  an 
electronically  addressed,  liquid  crystal  television  spatial  light 
modulator.  In  either  case,  the  data  was  coherently  illuminated  and 
Imaged  onto  the  first  HOE  (HI).  Fourier  transformation  with  a  lens 
then  produced  the  log  transformed  data  in  the  first  diffracted  order 
where  It  illuminated  the  second  HOE  (HI'  +  H2) .  The  second  HOE  and  a 
second  lens  then  performed  the  remainder  of  the  polar  formatting 
operation.  The  third  HOE  (H2')  corrected  the  output  phase,  and  a  final 
lens  was  used  to  perform  a  Fourier  transformation  and  produce  the 
desired  SAR  image. 
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V.  Experimental  Results 

The  multiple  HOE  system  was  used  experimentally  to  polar  format  and 
form  Images  from  SAR  data  of  multiple  point  scatterers.  Both  film  and 
a  liquid  crystal  television  (LCTV)  were  used  as  input  devices.  SAR 
data  was  written  on  the  input  devices  by  computer  generation  techniques 
and  the  output  image  was  magnified  and  detected  by  a  COHU  CCD  camera. 
The  overall  efficiency  of  the  system  was  about  1.6%;  however,  the 
output  Intensity  using  a  50  mW  laser  source  was  more  than  sufficient 
for  the  CCD  detector. 

Figure  7  shows  the  operation  of  the  system  at  each  stage  of  the 
formatting  and  processing  of  SAR  data  for  a  single  point  scatterer. 
For  this  experiment,  in  which  the  input  was  recorded  on  film  and  placed 
in  a  liquid  gate,  the  ability  of  the  multiple  HOE  system  to  properly 
format  the  data  and  provide  both  azimuth  and  range  compression  can 
clearly  be  seen. 

To  measure  system  performance  quantitatively,  SAR  data  representing 
eight  point  scatterers  In  a  single  line  was  generated  and  recorded  on 
film.  Because  the  data  was  recorded  on  a  bias  [as  in  Eq.  (24)],  an 
intense  DC  reference  point  and  high  sidelobes  along  the  x'  and  y'  axes 
are  present  in  the  output  Image  plane.  The  eight  point  test  image  was 
therefore  intentionally  placed  along  a  line  at  45°  to  these  axes  to 
reduce  the  effect  of  these  sidelobes.  In  Fig.  8(a),  the  image  Is  shown 
for  an  experiment  In  which  the  data  was  processed  without  polar 
formatting  (i.e.,  optical  Fourier  transformation  alone).  The  output  of 
the  multiple  HOE  system  is  shown  In  Fig.  8(b)  with  a  line  scan  through 
the  points  shown  in  Fig.  8(c).  Detector  noise  was  reduced  by  temporal 
averaging.  As  expected,  the  image  obtained  with  polar  formatting  shows 
far  superior  azimuth  compression. 
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Several  performance  characteristics  of  the  multiple  HOE  system  can 

be  determined  from  the  line  scan  in  Fig.  8(c).  First,  the  space- 

bandwidth  product  (SBWP)  can  be  measured.  Theoretically,  the  SBWP  is 

given  approximately  by  the  square  root  of  the  hologram  SBWP.6  For  the 
o 

(1  cm)  hologram  used  with  2048  by  2048  pixels,  a  SBWP  of  approximately 
45  by  45  would  be  predicted.  The  SAR  test  data  was  computed  such  that 
the  image  points  are  spaced  by  10  diffraction-limited  impulse  response 
widths.  The  system  modulation  transfer  function  (MTF),  specified  by 
the  envelope  of  the  line  scan,  falls  to  -3  dB  of  the  peak  value  at  62 
diffraction-limited  impulse  response  widths  along  the  diagonal  from  the 
DC  reference  point,  giving  an  experimental  SBWP  of  44  by  44. 

Second,  to  evaluate  the  system  phase  errors,  the  impulse  response 
widths  of  the  Images  of  the  eight  points  were  measured  and  found  to  be 
approximately  1.25  times  the  diffraction-limited  width  predicted  for  an 
ideal  system.  Not  all  of  this  increase  in  width  is  due  to  phase 
errors,  however.  Aperture  weighting  due  to  the  MTF  of  the  CGH  recorder 
is  partially  responsible  for  this  increase.  Another  part  of  the 
Increase  is  due  to  the  coordinate  transformation  itself.  This  aperture 
weighting  is  the  result  of  the  change  In  the  differential  element  of 
area  from  drd0  (rectilinear  r-0  coordinates)  to  rdrdfl  (polar 
coordinates)  and  can  also  be  predicted  by  evaluating  Eq.  (1)  using 
the  method  of  stationary  phase.11  Both  of  these  effects  can  be 
corrected  by  appropriate  Inverse  weighting  in  the  input  plane  of  the 
polar  formatting  system,  but  this  correction  was  not  implemented  in 
these  experiments. 

The  final  system  performance  characteristic  which  can  be  computed 
is  the  slgnal-to-nolse  ratio  (SNR).  The  system  SNR  Is  signal  dependent 
for  two  reasons.  First,  as  the  number  of  points  in  the  image 
increases,  the  available  system  dynamic  range  is  shared  by  more  points 
and  the  SNR  decreases.  Second,  SNR  variations  over  the  image  roughly 
follow  the  system  MTF  and  therefore  the  SNR  decreases  as  the  location 
of  the  point  moves  farther  from  the  DC  reference  point  (i.e.,  SNR 
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decreases  as  carrier  frequency  Increases).  For  the  point  of  the  eight 
point  image  which  was  nearest  to  the  DC  point,  a  SNR  of  34  dB  was 
measured. 

Real-time  operation  of  the  SAR  processor  was  demonstrated  by  using 

a  low  cost  Radio  Shack  Model  16-151  LCD  Pocketvision  TV  as  an  input 

spatial  light  modulator.  The  LCTV  was  modified  at  ERIM  for  this 

application  by  removing  the  polarizers  and  cementing  the  device  between 

19 

optical  flats  to  obtain  higher  optical  quality.  In  order  to  fully 
utilize  the  SBWP  of  the  LCTV,  the  input  imaging  optics  were  altered  to 
provide  3X  demagnification  of  the  LCTV  onto  the  first  HOE  HI.  The  LCTV 
was  driven  by  the  video  output  of  a  Data  Translation  DT2861  frame 
grabber  In  an  IBM  PC-AT  computer.  This  allowed  SAR  test  data  to  be 
placed  on  a  standard  video  signal  and  sent  to  the  LCTV.  A  polarizer 
oriented  at  90°  to  the  Input  laser  polarization  was  placed  just  before 
the  detector. 

The  output  of  this  real-time,  multiple  HOE,  SAR  polar  formatter  and 
processor  for  an  Input  representing  SAR  data  from  three  point 
scatterers  Is  shown  In  Fig.  9.  The  input  data  Is  shown  in  Fig.  9(a) 
and  the  output  image  In  Fig.  9(b).  Note  that  Fig.  9(b)  is  magnified 
with  respect  to  Fig.  8(b).  Because  the  SBWP  and  dynamic  range  of  the 
LCTV  are  less  than  that  of  the  film  input,  the  processor  performance  is 
reduced  somewhat.  The  smaller  SBWP  of  the  LCTV  can  be  seen  by 
comparing  the  MTF  of  the  real-time  processor  to  that  of  the  film-input 
processor  (see  Fig.  10).  The  real-time  processor  MTF  falls  off  much 
more  rapidly,  limiting  the  processor  SBWP  to  approximately  20  by  20. 
The  limited  dynamic  range  reduces  the  signal-to-blas  ratio  and  lowers 
the  SNR.  For  the  Image  of  the  test  point  nearest  the  DC  reference 
point,  a  SNR  of  approximately  18  dB  was  measured. 
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VI.  Conclusion 

A  multiple  HOE  system  was  designed,  constructed,  and  shown 
experimentally  to  perform  polar  formatting  for  the  real-time  optical 
processing  of  SAR  data.  An  electronically  addressed  SLM  was  used  for 
data  Input.  The  experimental  system  demonstrated  high  SNR,  low 
aberrations,  acceptable  efficiency,  and  a  SBWP  at  the  output  image  of 
44  by  44.  In  order  to  increase  the  processor  SBWP,  the  SBWP  of  the 
CGHs  used  must  be  increased  by  reducing  the  pixel  size  and  increasing 
the  CGH  size.  Using  a  (6.25  cm)2  CGH  with  a  1  /*m  spot  size,  a  250  x 
250  coordinate  transform  SBWP  could  be  achieved  which  would  also  match 
the  performance  of  current  real-time  input  and  output  devices. 

The  promising  results  obtained  in  this  research  provide  a  useful 
foundation  for  future  development  of  compact,  real-time  spotlight  SAR 
processors.  Further  work  in  this  area  should  first  be  aimed  at 
improving  the  system  performance  by  increasing  the  SBWP  and  reducing 
the  scattering  and  intermodulation  noise  in  the  COHOEs.  The  next  step 
would  be  the  design  of  a  more  compact  system  with  the  Fourier  transform 
lenses  incorporated  into  the  coordinate  transform  HOEs. 

This  research  was  supported  by  the  U.S.  Army  Research  Office. 
Portions  of  this  research  were  presented  at  the  Optical  Society  of 
America  Annual  Meeting  in  Seattle,  Washington,  October  20-24,  198613 
and  at  the  Society  of  Photo-Optical  Instrumentation  Engineers  Meeting 
on  Holographic  Optics:  Design  and  Applications  in  Los  Angeles, 
California,  January  13-14,  1988. 14 


28 


$>ERIM 


References 


1.  6.D.  Currie,  I.  Cindrich,  and  C.D.  Leonard,  "The  ERIM  TOPR  in 
Optical  Data  Processing,"  Proc.  SPIE  83,  8-13  (1976). 

2.  C.C.  Aleksoff,  I.J.  LaHale,  and  A.M.  Tai,  453  Optical -Hybrid 
Backprojectlon  Processing,"  Proc.  SPIE  422,  89-95  (1983). 

3.  M.  Haney  and  D.  Psaltls,  "Real-Time  Programmable  Acoustooptic 
Synthetic  Aperture  Radar  Processor,"  Appl.  Opt.  27,  1786-1796 
(1988). 

4.  0.  Bryngdahl,  "Optical  Map  Transformations,"  Opt.  Commun.  10,  164- 
168  (1974). 

5.  0.  Bryngdahl,  "Geometrical  Transformations  in  Optics,"  J.  Opt.  Soc. 
Am.  64.  1092-1099  (1974). 

6.  J.  Cederqulst  and  A.M.  Tai,  "Computer-Generated  Holograms  for 
Geometric  Transformations,"  Appl.  Opt.  23,  3099-3104  (1984). 

7.  W.J.  Hossack,  A.M.  Darling,  and  A.  Dahdouh,  "Coordinate 
Transformations  with  Multiple  Computer-Generated  Optical  Elements," 
J.  Modern  Optics  34.  1235-1250  (1987). 

8.  J.L.  Walker,  "Range-Doppler  Imaging  of  Rotating  Objects,"  IEEE 
Trans.  AES-16.  23-51  (1980). 

9.  M.S.  Welkowsky,  U.  Efron,  W.  Byles,  and  N.W.  Goodwin,  "Status  of 
the  Hughes  Charge-Coupled-Oevlce-Addressed  Liquid  Crystal  Light 
Valve,"  Opt.  Eng.  26,  414-417  (1987). 


29 


3>ERIM 


10.  S.H.  Lee,  S.C.  Esener,  M.A.  Title,  and  T.J.  Drabik,  "Two- 
Dimensional  Silicon/PLZT  Spatial  Light  Modulators:  Design 
Considerations  and  Technology,"  Opt.  Eng.  25,  250-260  (1986). 

11.  M.  Born  and  E.  Wolf,  Principles  of  Optics  (Pergamon,  New  York, 
1975),  pp.  752-754. 

12.  W.  Kaplan,  Advanced  Calculus  (Addison-Wesley,  Reading,  Mass., 
1959),  pp.  243-249. 

13.  J.N.  Cederquist  and  A.M.  Tai,  "Coordinate  Transformation  with 
Multiple  Computer-Originated  Holograms,"  J.  Opt.  Soc.  Am.  A  4, 
P121-P122  (1986). 

14.  J.N.  Cederquist,  A.M.  Tai,  and  M.T.  Eismann,  "Design  and 
Applications  of  Multiple  Hologram  Systems,"  Proc.  SPIE  883,  177-184 
(1988). 

15.  R.C.  Fairchild  and  J.R.  Fienup,  "Computer-Originated  Aspheric 
Holographic  Optical  Elements,"  Opt.  Eng.  21,  133-140  (1982). 

16.  H.  Bartelt  and  S.K.  Case,  "High  Efficiency  Hybrid  Computer 
Generated  Holograms,"  Appl.  Opt.  21,  2886-2890  (1982). 

17.  D.K.  Angell,  "Improved  Diffraction  Efficiency  of  Silver  Halide 
(Sensitized)  Gelatin,"  Appl.  Opt.  26,  4692-4702  (1987). 

18.  G.  Arfken,  Mathematical  Methods  for  Physicists  (Academic,  New  York, 
1966),  pp.  65-66. 

19.  A.M.  Tai,  "Low-cost  LCD  Spatial  Light  Modulator  with  High  Optical 
Quality,"  Appl.  Opt.  25,  1380-1382  (1986). 


30 


3>ERIM 


Data  from 
successive 
pulses  (0) 


Data  from  each  pulse  (r) 


i 


a)  rectilinear  format  b)  polar  format 


Fig.  1.  SAR  data  for  a  point  scatterer. 
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Fig.  2.  Optical  system  for  coordinate  transformation.  The  lens,  L, 

has  focal  length,  f,  and  Fourier  transforms  plane  P.  to  plane 

D  A 


32 


9/ERIM 


HI  lens  H1'andH2  lens  H2' 


Fig.  3.  Optical  system  for  polar  formatting. 
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.  COHOE  recording  method. 
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Input  HI  Output 


Fig.  6.  Experimental  multiple  HOE  system  for  polar  formatting  and  SAR 
image  formation. 
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Fig.  7.  Experimental  results  showing  polar  formatting  and  SAR  Image 
formation. 
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Fig.  8.  Experimental  results  for  SAR  data  for  eight  point  scatterers. 
2-D  Image  produced  (a)  without  polar  formatting  and  (b)  with 
polar  formatting,  (c)  Line  scan  through  point  images  in  (b) . 
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Fig.  9.  Experimental  results  with  LCTV  Input,  (a)  Data  for  three 
point  scatterers  which  was  input  to  LCTV.  (b)  Output  SAR 
Image. 
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Spatial  frequency  (Ip/mm) 


Fig.  10.  Relative  output  Image  intensity  versus  spatial  frequency  of 
Input  data  for  film  and  LCTV  input  devices. 
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APPENDIX  B 

ITERATIVE  DESIGN  OF  A  HOLOGRAPHIC  BEAMFORMER 
(accepted  by  Applied  Optics,  1989) 

Abstract 

A  design  method  is  presented  for  computing  the  phase  functions  of 
an  energy-efficient  system  using  two  holographic  elements  for 
converting  a  Gaussian  beam  into  a  uniform  beam  with  rectangular  support 
in  the  far-field  of  the  source.  The  method  Is  based  on  a  modification 
of  the  Gerchberg-Saxton  algorithm  which  Includes  an  x  -  y  separability 
constraint  on  the  phase  of  one  of  the  holographic  elements.  A 
beamforming  system  was  fabricated  using  this  method,  and  experimental 
results  were  obtained  which  support  the  design  approach. 

I.  Introduction 

In  applications  such  as  laser  radar  and  active  coherent  imaging,  it 
is  often  desired  to  uniformly  illuminate  a  rectangular  area  which  is  in 
the  far-field  of  the  source.  When  using  a  single  mode  laser  with  a 
Gaussian  (TEMQ0)  Intensity  distribution  as  a  source,  an  optical  device 
is  required  to  modify  the  near-field  amplitude  profile  of  the  beam  such 
that  its  phase  and  Intensity  become  uniform  upon  propagation  to  the 
far-field. 

Several  methods  have  been  investigated  for  beamforming  in  the  near- 

field.  The  most  straightforward  approach  Is  through  the  use  of 

apodlzation  and  truncation;  however,  this  approach  is  fundamentally 

energy  Inefficient.  For  energy  efficiency,  reflective1  or 
2  3 

refractive  '  surfaces  have  been  used  to  redistribute  the  energy  of  a 
Gaussian  beam  to  a  near-field  uniform  beam  of  circular  support. 
Unfortunately,  these  methods  generally  require  complicated  aspheric 
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surfaces  which  are  difficult  to  fabricate,  and  are  also  applicable  only 
In  a  limited  set  of  situations.  A  more  general  Implementation  approach 
Is  to  use  diffractive  optics  In  the  form  of  computer-generated 

4 

holograms.  The  conversion  of  a  Gaussian  beam  into  a  uniform  beam  of 
rectangular  support,  for  example,  has  been  demonstrated  using  this 

approach. ^ 

Diffractive  techniques  have  also  been  considered  to  obtain  a 
practical  beamformer  which  Is  useful  in  shaping  a  far-fleld  beam.  One 
system  uses  binary  phase,  Interlaced  gratings  to  provide  shaping  in  one 
dimension.6  In  this  system,  the  grating  pattern  provides  both 
intensity  weighting  and  the  required  phase  shifts  to  produce  a  1-D 
sinc(x)  =  (irx)-1  sin(irx)  distribution  In  the  near  field.  Upon 
propagation  to  the  far-field,  a  1-D  rect(x)  beam  profile  is  formed.  A 
similar  system  could  be  used  to  obtain  a  uniform  far-field  distribution 
with  a  circular  region  of  support  by  producing  a  2-D  J^xrJ/CTr) 
distribution  in  the  near-field.  The  system  is  fairly  efficient  in 
producing  a  JjUrJ/Crr)  distribution  even  though  weighting  is  used  to 
generate  the  desired  profile  due  to  the  similarity  between  the  Gaussian 
and  JjCirrJ/ffr)  functions.  System  efficiency,  however,  will  be  very 
poor  if  the  approach  Is  extended  to  producing  a  uniform  far-field  beam 
of  2-D  rectangular  support  since  the  necessary  near-field 
sinc(x)s1nc(y)  distribution  has  Its  energy  concentrated  near  the  axes. 
As  shown  In  the  Appendix,  the  efficiency  achievable  with  this 
apodization  approach  in  generating  a  sine  (x)sinc  (y)  intensity 
distribution  Is  only  5.8%.  This  efficiency  does  not  include  losses 
from  truncation  of  the  Gaussian  beam. 

To  obtain  greater  energy  efficiency,  we  investigated  the  use  of  a 
holographic  beamforming  system  to  provide  the  energy  redistribution 
Into  a  near-field  sine  (x)slnc  (y)  Intensity  profile  with  the  correct 
phase  without  the  use  of  intensity  weighting.  The  beamformer  concept 
Is  described  In  Section  II.  In  Section  III,  several  design  procedures 
are  described  for  synthesizing  the  required  phases  for  the  holographic 
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elements.  Experimental  procedures  and  results  are  then  discussed  in 
Section  IV.  Section  V  gives  conclusions  and  suggestions  for  further 
research. 

II.  Beamformer  Concept  Using  Two  Separated  Holographic  Elements 

A  single  holographic  optical  element  does  not  provide  enough 
degrees  of  freedom  to  produce  a  far-field  beam  with  both  uniform 
Intensity  and  phase  without  using  intensity  weighting.  In  order  to 
provide  both  the  energy  redistribution  and  phase  correction,  two 
separated  holographic  elements  are  required.  The  beamformer  concept 
therefore  uses  a  holographic  element  to  redistribute  the  energy  of  the 
Input  beam  into  the  necessary  near-field  Intensity  profile,  and  then 
uses  a  second  holographic  element  to  place  the  correct  phase  onto  this 
beam.  By  relying  solely  on  energy  redistribution  and  phase  correction, 
the  system  Is  limited  In  efficiency  only  by  the  maximum  achievable 
diffraction  efficiency  of  the  holographic  elements,  which  can  be  very 
high. 

The  basic  beamformer  system  Is  Illustrated  in  Fig.  1.  The 
truncated  Gaussian  beam  from  a  single  transverse  mode  laser  is  incident 
on  a  holographic  element  (HI)  which  redistributes  the  energy  into  a 
desired  near-field  intensity  distribution.  Then,  a  second  holographic 
element  (H2)  places  the  desired  near-field  phase  distribution  onto  the 
beam  such  that,  upon  propagation  to  the  far-field,  a  uniform  Intensity 
and  phase  distribution  of  approximately  rectangular  support  results. 
After  determining  a  near-field  amplitude  distribution  which  yields  the 
desired  far-field  characteristics,  the  design  of  the  beamforming  system 
consists  of  computing  the  required  phases  for  HI  and  H2  to  produce  the 
desired  near-field  beam. 

Since  the  illuminating  beam  Is  coherent,  a  Fourier  transform 
relationship  exists  between  the  near-field  and  far-field  amplitude 
distributions.7  To  produce  a  uniform  far-field  distribution  of 
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rectangular  support,  denoted  as  rect(x)  rect(y),  a  near-field  sinc(x) 
sinc(y)  distribution  Is  required.  This  distribution,  however,  is  of 
infinite  extent,  and  Is  therefore  unrealizable.  To  obtain  a  realizable 
near-fleld  solution,  the  Ideal  solution  must  be  windowed  by  some 
function  which  is  zero  outside  the  area  of  the  beamformer  elements. 

In  choosing  the  optimum  window  function,  one  inevitably  encounters 
a  tradeoff  between  uniformity  and  the  ability  to  obtain  sharp  edges  in 
the  far-field  due  to  the  finite  aperture  size  available.  A  rectangular 
window  (i.e.,  truncation)  In  the  near-field  results  in  the  sharpest 
far-field  edges,  but  with  very  significant  ringing.  In  order  to  reduce 
the  ringing,  the  window  function  must  roll  off  more  uniformly  to  zero, 
but,  in  doing  this,  the  far-field  edges  become  less  sharp.  For  the 
purpose  of  this  investigation,  a  linear-weighted  window  with  a  cutoff 
at  the  sixth  zero  of  the  sine  function  was  arbitrarily  chosen.  The 
desired  near-field  amplitude  distribution,  AN(x,y),  is  given  by 

AN(x,y)  =  sinc(6x/x0,  6y/yJ  (l  -  lxl/xj(l  -  lyl/yj  ,  (1) 

where  the  second  holographic  element  (H2)  size  is  2xQ  by  2yQ.  The 
truncated  Gaussian  Input  beam  was  chosen  to  have  a  10%  truncation  loss. 
The  modulus  of  this  beam  is  shown  in  Fig.  2(a).  The  magnitude  of  the 
distribution  from  Eq.  (1)  Is  shown  in  Fig.  2(b),  and  its  corresponding 
far-field  beam  modulus  In  Fig.  2(c). 

III.  Holographic  Element  Design 

Now  that  the  desired  near-fleld  distribution  has  been  chosen,  a 
holographic  system  must  be  designed  which  can  produce  this  distribution 
from  a  uniform  phase,  truncated  Gaussian  input.  Referring  to  Fig.  1, 
holographic  element  HI  provides  the  energy  redistribution  from  the 

O 

truncated  Gaussian  beam  Into  a  beam  with  a  weighted  and  truncated  sine 
Intensity  profile.  In  this  section,  several  methods  are  discussed  for 
determining  the  phase  of  this  hologram.  The  second  holographic  element 
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H2  corrects  the  phase  of  the  output  beam  of  HI,  and  is  discussed  at  the 
end  of  this  section. 

A.  Design  Based  on  Stationary  Phase  Approximation 

The  first  design  method  to  be  discussed  Is  based  on  the  method  of 
8  9 

Bryngdahl  '  for  producing  geometric  transformations  with  holographic 
elements.  The  basic  geometry  is  shown  in  Fig.  3(a),  where  the  input 
field  at  Pj  has  an  amplitude  distribution  a^x.y)  with  uniform  phase. 
It  is  desired  to  add  a  holographic  element  with  phase  0(x,y)  in  plane  Px 
such  that  a  given  output  amplitude  distribution  a2(u,v)  results  in  the 
back  focal  plane  P2  of  lens  L.  The  approximate  relationship  between  the 
two  planes  is  given  by  the  Fourier  transform: 

a2(u,v)  =  ’  if  f  f  ai(x»y)  exp[i*(x,y)  -  ^  (xu  +  yv)]  dxdy  (2) 

where  X  is  the  wavelength,  f  is  the  focal  length  of  lens  L,  and  (x,y) 
and  (u, v)  are  the  coordinates  in  planes  Pj  and  P2,  respectively.  To 
simplify  this  system,  note  that  by  placing  the  input  plane  Pj  in  contact 
with  the  lens  l,  as  in  Fig.  3(b),  the  output  distribution  only  changes 
by  a  quadratic  phase: 

a2(u,v)  «  exp[  ^  [u2  +  v2]]  a2(u,v).  (3) 

Further,  as  shown  in  Fig.  3(c),  the  quadratic  phase  of  the  lens  may  be 
Included  In  the  holographic  element  by  defining  a  new  phase  where 

#’(x,y)  3  #(x,y)  -  p  (x2  +  y2]  .  (4) 

Therefore,  the  necessary  phase  of  the  first  beamforming  hologram  (HI)  in 
Fig.  1  can  be  derived  from  ^(x,y). 
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By  approximating  the  Fourier  transform  Integral  In  Eq.  (2)  by  the 
method  of  stationary  phase,10  It  can  be  shown  that  a  point  (x,y)  in  the 
object  plane  Is  mapped  to  a  point  (u,v)  In  the  Fourier  plane  given  by 
the  relation: 


If  the  desired  mapping  is  defined  by 


the  necessary  phase  to  produce  this  mapping  is  the  solution  of  the 
partial  differential  equations: 


lx  *  Xf  g(x'y)' 
ly  =  X?  h(x’y)  * 


(7) 


For  conversion  of  a  truncated  Gaussian  beam  into  a  weighted  sine 
distribution,  both  distributions  are  separable  functions.  Because  of 
this,  the  partial  differential  equations  in  Eq.  (7)  become  ordinary 
differential  equations  In  x  and  y.  If  the  mapping  functions  g(x)  and 
h(y)  are  bounded,  the  required  coordinate  transforms  can  be  found 
independently  in  the  x  and  y  directions.  Once  the  mapping  functions  are 
determined,  the  phase  of  the  hologram  is  given  by 


v  y 

J  g(x ' )dx '  +  j  h (y ' ) dy '  -  ff  (*2  +  y2) • 


(8) 


Note  that  the  quadratic  phase  of  the  lens  has  been  added. 
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The  method  of  conservation  of  energy11  can  be  used  In  both  the  x  and 
y  directions  to  determine  the  mapping  functions.  Let  one  dimension  of 
the  Incident  Gaussian  Intensity  be  represented  by 


Vx) 


'I  exp(-x2/2a2)  Ixl  i  xQ 
.0  Ixl  >  x0, 


(9) 


where  I  Is  the  peak  Intensity,  xQ  is  the  half-width  of  the  truncated 
beam,  and  a  Is  a  constant.  From  Eq.  (1),  the  desired  near-field  1-D 
Intensity  distribution  is 


Vu> 


ho  s1"cZ(6u/uo)(1  -  lul/uo)2 
.  0 


lul  i  u0 
lul  >  uo  . 


(10) 


The  mapping  can  be  determined  with  the  aid  of  Fig.  4.  All  of  the 

energy  bounded  by  0  i  x*  S  x  from  the  Gaussian  intensity  profile  is 

2 

mapped  into  a  region  bounded  by  0  $  uM  u  in  the  sine  intensity 
distribution.  This  can  be  described  mathematically  by 


x  u 

lo  J  exP(~x,2/2a2)  dx1  =  INo  J  s1nc2(6u7uQ)  (1  -  I u '  I /uQ)2  du'  (11) 
o  o 


for  Ixl  <  xQ  and  lul  <  uQ.  The  peak  intensity  in  the  two  planes  is 
related  by  conservation  of  the  total  energy  in  both  distributions,  or 


J  Ig(x)  dx  =  J  IN(u)  du  .  (12) 
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Inserting  the  functions  given  In  Eqs.  (9)  and  (10)  and  solving  for  INo 
gives 


0 

J  exp(-x2/2 a2)  dx 


which  can  be  evaluated  numerically. 

The  mapping  function  g(x)  can  be  determined  by  numerically  solving 
Eq.  (11)  to  determine  a  set  of  points  u  as  a  function  of  x.  This  is 
Implemented  by  computing  the  left-hand-side  of  Eq.  (11)  for  a  given 
value  of  x,  and  then  computing  the  right-hand-side  for  Increasing  values 
of  u  until  the  equality  Is  satisfied.  Due  to  the  separability  of  the 
Gaussian  and  sine  functions,  h(y)  can  be  similarly  obtained.  With  the 
mapping  functions  defined,  Eq.  (8)  can  then  be  numerically  Integrated  to 
determine  the  needed  phase  function  for  the  light  redistribution 
hologram. 

For  this  Investigation,  a  symmetric  Gaussian  beam  truncated  to 

p 

contain  90%  of  its  energy  within  a  (10.24mm)  area  was  used  as  an  input. 

2 

The  beamformer  was  designed  to  produce  the  weighted  sine  intensity  of 
Eq.  (10)  with  xQ  and  yQ  equal  to  5.12mm.  The  separation  between  the 
holographic  elements  was  375mm.  After  evaluating  Eq.  (13)  numerically 
for  INq,  the  mapping  function  was  determined  by  balancing  the  integrals 
of  Eq.  (11)  for  100  (x,u)  pairs.  The  resulting  function  was 

Interpolated  to  1024  points  using  a  cubic  spline  Interpolation  and 
numerically  Integrated  using  Eq.  (8)  to  compute  the  required  holographic 
element  phase. 
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The  stationary  phase  design  method  was  tested  by  computer 
simulation.  This  was  done  by  creating  a  256  x  256  array  with  Gaussian 
Intensity  and  the  numerically  determined  phase  (without  the  quadratic 
phase),  padding  It  with  zeros  to  a  512  x  512  size,  and  taking  a  2D  FFT 
to  give  the  near-field  distribution.  The  simulated  far-field 
distribution  was  computed  by  replacing  the  simulated  near-field  phase 
with  the  Ideal  phase  (alternating  r  phase  shifts  of  the  sine  function), 
and  taking  a  2D  FFT.  The  resulting  near-field  and  simulated  far-field 
moduli  are  shown  in  Fig.  5.  It  can  be  seen  by  comparing  to  Fig.  2  that, 
while  the  stationary  phase  approximation  was  sufficient  to  provide  the 
general  shape  of  the  near-field  distribution,  it  was  clearly 
insufficient  for  accurately  producing  the  zeros  of  the  distribution,  and 
results  In  unacceptable  ringing  and  high  sldelobes  in  the  far-field. 

B.  Design  Using  Iterative  Fourier  Transform  Algorithms 

To  Improve  the  design  of  the  previous  section,  a  second  design 

method  was  Investigated.  Because  of  the  Fourier  transform  relationship 

between  the  two  planes  In  Fig.  3(a),  the  design  problem  Is  equivalent  to 

that  of  computing  an  object  phase  (CGH  phase  not  Including  qu^Jratic 

phase)  for  a  known  object  modulus  (truncated  Gaussian)  and  a  known 

Fourier  modulus  (weighted  sine).  This  synthesis  problem  Is  equivalent  to 

12 

that  of  phase  retrieval  from  a  pair  of  intensity  measurements,  in 

which  the  complex  fields  In  two  domains  related  by  a  Fourier  transform 

are  reconstructed  from  Incomplete  information  in  each  domain.  This 

problem  has  gained  much  attention  In  the  past,  in  applications  such  as 

electron  microscopy,  and  Is  typically  solved  using  an  iterative  Fourier 

transform  algorithm.  For  this  specific  case,  in  which  the  object  and 

13 

Fourier  moduli  are  known,  the  Gerchberg-Saxton  algorithm  allows 
reconstruction  of  the  phase. 

A  block  diagram  of  the  Gerchberg-Saxton  algorithm  is  shown  in  Fig. 
6.  The  algorithm  consists  of  an  iterative  loop  having  four  steps:  (1) 
the  known  object  modulus  with  an  Initial  phase  estimate  is 


49 


22ERIM 


transformed  to  the  Fourier  domain,  (2)  the  Fourier  modulus  IGI  Is 
replaced  with  the  known  or  desired  modulus,  (3)  the  resulting  Fourier 
data  Is  transformed  back  to  the  Image  domain,  and  (4)  the  object  modulus 
Ig'l  Is  replaced  with  the  known  object  modulus.  The  resulting  object 
data  Is  then  input  to  step  (1),  and  the  cycle  Is  repeated  until 
convergence  Is  achieved.  This  algorithm  has  been  used  to  design 
holographic  beamformers14, 15  and  computer-generated  Fourier 
holograms.16-18 

A  2-D  Iterative  synthesis  of  the  first  holographic  element  phase 
using  this  algorithm  was  attempted.  The  results  were  found  to  be  very 
sensitive  to  the  Initial  phase  estimate  used  to  start  the  algorithm.  It 
was  found  that  by  starting  with  a  random  phase,  the  resulting 
redistributed  modulus  Is  approximately  a  weighted  sine,  but  with 
randomly  distributed  error,  as  shown  in  Fig.  7.  The  resulting  CGH 
phase,  though  an  approximate  solution,  has  significant  random  features 
which  vary  for  each  different  random  input.  This  lack  of  repeatability 
was  considered  unacceptable  for  a  synthesis  procedure. 

Better  results  were  obtained  by  using  the  stationary  phase  solution 
as  an  Initial  phase  estimate.  After  300  iterations,  the  algorithm 
converged  to  the  solution  whose  simulated  near-field  and  far-field 
moduli  are  shown  in  Fig.  8.  Although  the  algorithm  was  Initialized  with 
a  separable  object  modulus,  object  phase,  and  Fourier  modulus,  the 
resulting  object  phase  became  less  separable  as  the  number  of  iterations 
Increased.  This  Is  attributed  to  finite  word  length  effects  in  the 
digital  computation. 

This  lack  of  separability  In  the  phase  solution  Is  an  undesirable 
feature  of  the  Gerchberg-Saxton  result.  The  solution  would  not  be 
expected  to  be  optimal,  since  the  coordinate  transform  problem  Is 
exactly  separable.  Also,  a  separable  solution  would  be  easier  to  encode 
as  a  CGH,  because  It  would  make  the  Interpolation  to  the  2048  x  2048 
holographic  element  size  less  computationally  intensive.  For  these 
reasons,  It  was  desired  to  find  a  separable  phase  solution. 
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The  Iterative  algorithm  was  therefore  modified  to  enforce  a 
separability  constraint  on  the  object  phase.  This  hybrid  algorithm  is 
shown  in  Fig.  9.  An  initial  phase  estimate  is  used  to  start  the 
Gerchberg-Saxton  algorithm,  and  a  specified  number  of  iterations  are 
performed.  Then  a  1-D  slice  of  the  current  phase  estimate  is  extracted, 
and  an  outer  product  is  performed  to  produce  a  2-D  separable  phase 
estimate.  This  is  used  to  start  another  set  of  iterations,  and  the  loop 
is  repeated  until  an  acceptable  solution  is  found  by  examining  the 
simulated  far-fleld  intensity.  After  each  separability  constraint  is 
imposed,  the  phase  estimate  improves,  so  that  the  iterative  algorithm 
converges  with  fewer  iterations.  Each  successive  phase  solution  is 
closer  to  being  separable,  so  that  overall  convergence  to  a  good 
separable  phase  solution  would  be  expected. 

A  phase  solution  was  generated  by  performing  5  loops  of  the  hybrid 
algorithm  with  a  successively  decreasing  number  of  Gerchberg-Saxton 
iterations  per  hybrid  loop  (100,  50,  25,  10,  and  5).  The  results  are 
shown  in  Fig.  10.  When  compared  to  the  nonseparable  phase  solution 
generated  by  the  Gerchberg-Saxton  algorithm  alone  [Fig.  9(a)],  the 
hybrid  solution  has  slightly  better  uniformity  and  lower  sidelobes  in 
the  simulated  far-fleld  beam  for  the  same  total  number  of  iterations. 
Since  it  is  separable,  it  is  also  easier  to  encode  as  a  CGH. 

C.  Design  of  Second  Holographic  Element 

To  complete  the  beamformer  system  design,  the  phase  of  the  second 
holographic  element  H2  must  be  specified.  The  purpose  of  H2  is  to  place 
on  the  beam  the  desired  near-field  phase  which  Is  a  simple,  known 
function  consisting  of  alternating  t  phase  shifts  at  the  zeroes  of  the 
slnc(x)  slnc(y)  distribution.  H2  must  also  compensate  for  the  non- 
uniform  phase  of  the  beam  arriving  from  HI.  It  is  therefore  necessary 
to  determine  this  phase  function,  which  can  be  done  by  numerically 
propagating  the  complex  amplitude  at  the  plane  of  HI  to  that  of  H2,  and 
extracting  the  phase.  Since  this  may  require  some  approximations  and 
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would  Ignore  any  aberrations  In  the  system,  it  may  be  limited  in  its 
accuracy.  For  this  reason,  It  may  be  necessary  to  measure  the  phase  at 
the  plane  of  H2  with  an  interferometric  system.  Once  this  Is  done,  the 
required  phase  of  H2  can  also  be  encoded  in  a  C6H.  These  two  methods 
may  be  very  difficult  to  perform  experimentally.  First,  it  may  be 
difficult  to  compute  or  measure  the  phase  at  the  plane  of  H2  within  the 
needed  accuracy  since  the  beam  phase  changes  very  rapidly  in  this 
region.  Second,  even  if  the  phase  can  be  determined,  the  tolerances  for 
aligning  H2  to  HI  would  be  very  tight.  Small  errors  in  alignment  or 
phase  compensation  would  lead  to  unacceptable  results  in  the 
beamforming. 

A  third  method  for  producing  H2  Is  by  interfering  the  output  of  HI 
with  a  beam  carrying  the  desired  near-field  phase.  Upon  reconstruction, 
the  desired  phase  Is  imposed  on  the  beam  after  it  is  diffracted  by  H2. 
Ideally,  the  diffraction  efficiency  across  the  entire  phase  correction 
hologram  should  be  constant  in  order  for  the  diffracted  output  to 
maintain  the  desired  profile.  To  achieve  uniform  diffraction  efficiency 
with  a  recording  beam  (coming  from  HI)  having  a  nonuniform  intensity 
distribution  requires  a  highly  non-linear  recording  process. 
Unfortunately,  such  a  non-linear  recording  will  result  in  high 
Intermodulation  noise.  To  avoid  the  problem  with  Intermodulation  noise, 
the  recording  has  to  be  kept  linear.  A  linear  recording,  however, 
produces  a  diffraction  efficiency  that  Is  proportional  to  the  input 
intensity  distribution.  The  result  is  a  squaring  of  the  input  intensity 
profile.  This  squaring  effect  can  be  compensated  by  designing  the  light 
redistribution  hologram  to  produce  an  Intensity  distribution  equal  to 
the  square  root  of  the  desired  intensity  distribution.  Another  drawback 
of  restricting  the  recording  to  the  linear  region  Is  the  low  average 
diffraction  efficiency  that  can  be  obtained. 
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IV.  Experiment 

To  verify  the  results  of  the  computer  simulation,  HI  was  produced 
using  the  phases  synthesized  by  both  the  stationary  phase  and  the  hybrid 
algorithm  design  methods.  The  conclusions  which  were  drawn  earlier  by 
simulation  of  the  stationary  phase  design  method  are  supported  by 
experimental  results.  A  CGH  was  produced  by  adding  the  phase  solution 
for  HI  from  this  method  to  an  appropriate  quadratic  phase,  encoding  the 
phase  as  a  real,  non-negative  amplitude  transmittance  function  on  a 
carrier,  and  recording  the  distribution  on  film  using  a  rotating  drum 
laser  beam  recorder.4  A  hologram  size  of  (1.024  cm)^  was  used  with  a 
5pm  pixel  size  to  produce  a  hologram  space-bandwidth  product  of  2048  x 
2048.  The  hologram  was  mounted  between  1/4  inch  microflats  using  UV 
curing  cement  In  order  to  eliminate  phase  errors  due  to  thickness 
variations  in  the  film.  When  illuminated  with  the  Gaussian  beam,  the 
redistributed  near-field  Intensity  profile  occurs  in  the  first 
diffracted  order  in  a  plane  375mm  away,  and  is  shown  in  Fig.  11. 
Because  the  diffraction  efficiency  of  the  CGH  was  very  low  (~  1%),  the 
intensity  distribution  was  measured  with  a  scanning  fiber  optic  probe 
which  was  coupled  Into  a  photomultiplier  tube  detector.  The  slight 
asymmetry  in  the  distribution  is  due  to  the  MTF  of  the  laser  writer 
( 1 . e . ,  higher  spatial  frequencies  are  diffracted  less  strongly). 

Significantly  better  results  were  obtained  with  the  hybrid  algorithm 
design  method.  A  CGH  with  the  designed  phase  function  was  produced  as 
above  using  the  laser  writer.  The  intensity  distribution  of  the 
diffracted  Gaussian  beam  at  plane  H2  is  shown  in  Fig.  12(a).  Since  the 
result  of  the  hybrid  algorithm  design  method  is  clearly  better  than  that 
of  the  stationary  phase  method,  the  CGH  designed  using  the  former 
approach  was  used  for  HI.  The  CGH  recorded  on  film  as  an  absorption 
hologram,  however,  Is  too  low  In  diffraction  efficiency  (~  1%)  to  be 
practical.  In  order  to  Improve  the  diffraction  efficiency,  the  CGH  was 
copied  onto  a  volume  phase  material  with  a  high  carrier  frequency,  which 
then  becomes  known  as  a  computer-originated  holographic  optical  element 
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(COHOE).  The  COHOE  was  recorded  as  shown  in  Fig.  13  by  exposing  a 

Kodak  649F  spectroscopic  plate  using  a  HeNe  laser  (X  =  0.6328  urn),  and 

processing  using  the  silver  halide  (sensitized)  gelatin  process.  A 

diffraction  efficiency  of  approximately  43%  was  achieved,  v'hich  could  be 

21 

increased  by  using  the  dichromated  gelatin  process  end  an  argon  laser. 
When  reconstructed  in  the  same  manner  as  the  CGH,  the  resulting  near- 
field  Intensity  of  the  C0H0E  is  shown  in  Fig.  12(b). 

A  note  should  be  made  here  about  the  methods  used  to  obtain  the 
intensity  line  scans.  The  Intensity  distributions  of  the  output  of  the 
CGHs  [Figures  11(b)  and  12(e)]  were  obtained  with  the  scanning  fiber 
optic  probe  since  the  diffraction  efficiency  was  so  low.  For  the  output 
of  the  COHOE  [Figure  12(f)],  however,  more  light  was  available  for 
detection,  so  the  line  scans  were  obtained  by  using  a  CCD  camera  and 
grabbing  video  lines  with  a  digital  oscilloscope.  The  differing  methods 
of  detection  account  for  the  different  characteristics  of  the  line 
scans.  In  particular,  tue  pixel  size  of  the  fiber  optic  probe  is  larger 
than  that  of  the  CCD  camera,  resulting  In  more  spatial  integration  which 
smooths  the  output  of  the  probe. 

The  second  hologram  H2  was  produced  using  the  third  method  described 
at  the  end  of  Section  III  by  interfering  the  output  of  the  first 
hologram  HI  with  a  beam  carrying  the  desired  near-field  phase  of  the 
sine  distribution.  This  desired  phase  (alternating  r  phase  shifts)  was 
encoded  and  recorded  as  a  CGH  on  film.  This  step  was  very  sensitive  to 
alignment  errors,  since  phase  compensation  has  to  occur  at  the  right 
plane  and  the  desired  r  phase  shifts  must  line  up  almost  exactly  with 
the  zeroes  of  the  sine  distribution.  Very  slight  misalignment  errors 
were  found  to  result  In  poor  beamforming. 

An  added  problem  Involves  the  linearity  of  the  recording.  Since  the 
sldelobes  of  the  near-field  distribution  are  much  lower  than  the 
malnlobe,  a  linear  recording  over  a  very  wide  dynamic  range  is  required. 
Relatively  small  nonl inearities  result  in  Incorrect  sidelobe-to-mainlobe 
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Intensity  ratios,  which  correspond  to  either  excessive  ringing  or  edge 
blurring  In  the  far-fleld  beam.  Keeping  the  recording  very  linear 
severely  limits  the  overall  diffraction  efficiency  that  can  be  obtained. 
The  resulting  holographic  element  In  our  experiment  had  a  diffraction 
efficiency  of  only  5.2%. 

The  results  from  the  complete  beamforming  system  are  shown  in  Fig. 
12(c)  and  (d) .  In  the  near-field,  the  squaring  effect  of  the  second 
hologram  can  be  readily  seen  by  the  relative  sidelobe  intensities.  In 
the  far-fleld,  the  beam  Is  relatively  uniform  with  very  low  sidelobes 
and  sharp  edges,  although  there  Is  still  some  ringing  In  the  edges  and 
corners.  The  entire  distribution  also  seems  to  be  weighted  slightly 
higher  In  one  corner.  This  Is  due  to  the  presence  of  nonlinearities  in 
the  recording  and  small  misalignment  errors.  The  fringes  seen  in  the 
far-field  beam  are  due  to  Interference  in  the  faceplate  of  the  CCD 
camera,  and  are  not  present  In  the  actual  output  beam  of  the  beamforming 
system.  The  lira  scans  [Figures  12(g)  and  12(h)]  were  obtained  from  the 
video  output  of  the  CCD  camera. 

V.  Discussion 

In  order  to  obtain  a  highly  efficient,  uniform  far-field  beamforming 
device,  the  system  must  rely  solely  on  redistribution  of  energy. 
Several  methods  were  investigated  for  computing  the  necessary  phase 
functions  of  a  holographic  system  for  this  application.  It  was  found 
that  a  method  based  on  using  a  stationary  phase  approximation  to  the 
Fourier  transform  Integral  was  inadequate  for  producing  good  zeroes  in 
the  near-field  distribution.  Significantly  better  results,  however, 
were  achieved  using  a  hybrid  algorithm  which  combined  an  x-y 
separability  constraint  with  an  Iterative  Fourier  transform  algorithm. 
Further  Improvements  may  also  be  made  by  using  alternate  algorithms  such 
as  the  Input-output  algorithm  that  speeds  convergence.  '  Another 
possibility  would  be  to  use  a  two-step  design  procedure  which  would 
optimize  the  beam  In  the  far-fleld  plane  rather  than  specifying  the 
desired  beam  at  H2. 
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Experimental  results  demonstrated  the  ability  of  the  holographic 

system  to  provide  accurate  beamshaping.  To  obtain  uniformly  high 

efficiency  across  the  second  holographic  element,  it  must  be  fabricated 

as  a  COHOE  by  computing  or  measuring  the  necessary  phase  as  described  in 

Section  III.  While  this  puts  very  high  demands  both  on  alignment  and 

the  ability  to  produce  precise  phase  compensation,  it  is  necessary  in 

order  to  make  an  efficient  far-field  beamforming  device.  By  doing  this, 

the  overall  efficiency  becomes  limited  only  by  that  of  the  holographic 

optical  elements  which  would  be  expected  to  be  >80%  with  existing 

materials  (i.e.,  dichromated  gelatin).  For  comparison,  using  intensity 

o 

weighting  to  obtain  the  near-field  weighted  sine  intensity  distribution 
of  Eq.  (10),  the  maximum  efficiency  that  can  be  achieved  is  only  8%  (see 
Addendum).  Neither  of  these  efficiencies  include  losses  due  to 
truncation  of  the  input  Gaussian  beam. 

This  research  was  sponsored  by  the  U.S.  Army  Research  Office.  The 
authors  would  like  to  thank  J.R.  Fienup  for  helpful  advice  on  the  use  of 
the  Gerchberg-Saxton  algorithm. 
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ADDENDUM 

o 

Since  both  the  truncated  Gaussian  and  weighted  sine  intensity 

distributions  are  separable  functions,  the  efficiency  In  apodizing  a  2-D 

2 

Gaussian  Into  a  weighted  sine  beam  profile  Is  the  square  of  the  1-D 
efficiency.  Let  the  1-D  truncated  Gaussian  intensity  be  described  by 


IG(X) 


[i0  exp(-x2/2a2) 


and  the  1-D  weighted  sine 


1  0 

2 

Intensity  by 


Ixl  S  xQ 
Ixl  >  x0 


!n(x) 


INq  s1nc2(Nx/xQ)(l  -  lxl/xQ)2 
.  0 


1  X I  S  x0 
Ixl  >  x0 


(Al) 


(A2) 


where  2xq  is  the  window  width  and  (N-l)  Is  the  number  of  sidelobes  on 
each  side  of  the  sine  function  inside  the  window. 


For  efficient  apodlzatlon,  the  intensities  are  equalized  at  x  =  0: 


I 


No  = 


(A3) 


In  order  to  exactly  produce  the  distribution  In  Eq.  (A2)  without  energy 
redistribution,  the  a  parameter  of  the  Gaussian  must  be  chosen  such  that 

IG(x)  l  IN(x)  for  Ixl  S  xQ.  (A4) 

The  greatest  efficiency  Is  obtained  when  the  Intensities  are  equal  at 
the  peak  of  the  (N-l)th  sldelobe,  which  Is  located  approximately  at 
x  =  [(N  -  1/2) /N]  xQ.  The  condition  Is  therefore  given  by 
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- 


l0  exp 


fW  -  1/21 2  Xo 
1  N  J  2*2J 


=  I. 


s1nc2[(N  -  1/2)] [l  -  [N-'n1/2]]  .  (A5) 


This  can  be  solved  for  the  required  a  parameter: 


a  = 


21n[4N2[N  -  1/2)  V] 


1/2 


(A6) 


The  1-D  apodlzatlon  efficiency,  7,  Is  equal  to  the  ratio  of  the 
energy  in  the  apodlzed  beam  to  that  in  the  incident  beam.  For  this 
case, 


|  0  IQ  sinc2[Nx/xQ]  [l  -  Ixl/xJ2  dx 

r,  =  — ^ - - -  (A7) 

j  0  IQ  exp£-x2/2a2]  dx 
"xo 

o 

where  a  is  given  by  Eq.  (A6) .  The  2-D  efficiency  is  17  .  For  the  case 

where  N  =  6,  the  Integrals  In  Eq.  (A7)  can  be  computed  numerically  to 
2 

obtain  rj  =  8.1%.  This  efficiency  can  be  improved  by  using  a  smaller  a 
and  clipping  the  higher  sldelobes  of  the  sine  intensity,  but  this  will 
result  In  some  degradation  In  the  far-field  beam.6  We  have  also 
calculated  the  efficiency  of  converting  a  Gaussian  beam  into  a  non- 
weighted  sine  intensity  distribution.  This  was  done  by  eliminating  the 
weighting  term  in  Eq.  (A2)  and  computing  the  efficiency  in  the  same 
manner  as  shown  above.  In  this  case,  the  efficiency  is  only  5.8%. 
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INPUT 

AMPLITUDE 


FAR-FIELD 

AMPLITUDE 


Fig.  1.  Holographic  beamformer.  Holographic  element  HI  redistributes 
the  energy  of  the  incident  beam  into  the  desired  near-field 
intensity  distribution.  Holographic  element  H2  places  the 
correct  phase  on  the  beam.  Upon  propagation  to  the  far-field, 
a  uniform  intensity  and  phase  results. 
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A 

V 


/ 


/ 


(a) 


a1(x,y)el<{>(x,y)  a2’(u,v) 


/ 


(b) 


a^x  ,y  )  e1^  {X,y)  a2’(u,v) 


/ 


(C) 


Fig.  3.  Equivalent  systems  for  holographic  redistribution  of  energy. 

In  (a),  the  complex  amplitude  distributions  are  related  by  a 
Fourier  transform.  If  the  holographic  element  Is  placed  in 
contact  with  the  lens,  as  In  (b),  the  Intensity  distribution 
at  the  output  plane  Is  unchanged.  In  (c),  the  quadratic  phase 
of  the  lens  Is  added  to  the  holographic  element. 
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Fig.  5.  Simulated  results  for  beamformer  designed  by  stationary  phase 
method,  (a)  Near-field  modulus;  (b)  far-fleld  modulus. 
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Fig.  6.  The  Gerchberg-Saxton  algorithm. 


07 


88-11487 


Fig.  7.  Near-field  modulus  obtained  using  Gerchberg-Saxton  algorithm 
with  a  random  initial  phase  estimate. 


88-11488 


(a) 


(b) 


Fig.  8.  Synthesis  results  using  Gerchberg-Saxton  algorithm  with 
initial  phase  estimate  from  stationary  phase  method,  (a) 
Near-field  modulus;  (b)  far-field  modulus. 
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SEPARABILITY  GERCHBERG  -  SAXTON 

CONSTRAINT  ALGORITHM 


Fig.  9.  The  hybrid  Iterative  algorithm.  The  separability  constraint 
Is  applied  after  each  set  of  Gerchberg-Saxton  iterations. 
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Fig.  10.  Synthesis  results  for  beamformer  designed  using  hybrid 

algorithm,  (a)  Near-field  modulus;  (b)  far-field  modulus. 
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(b) 


Fig.  11.  Experimental  results  for  CGH  computed  by  the  stationary  phase 
method,  (a)  Intensity  distribution  at  a  plane  375mm  from  the 
CGH;  (b)  line  scan  through  (a). 
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(a)  (b)  (c)  (d) 


Fig.  12.  Beamformer  experimental  results,  (a)  Intensity  distribution 
at  a  plane  375  mm  from  CGH;  (b)  same  result  from  COHOE;  (c) 
near-field  intensity  distribution  directly  after  H2;  (d)  far- 
fleld  Intensity  distribution;  (e)-(h)  are  line  scans  through 
the  distributions  in  (a)-(d),  respectively. 
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COLLIMATOR 


Fig.  13.  COHOE  fabrication.  The  desired  first  diffraction  order  of  the 
CGH  is  Imaged  onto  the  COHOE  plane  and  interfered  with  a 
reference  beam. 
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COORDINATE  TRANSFORMATIONS  REALIZABLE  WITH 
MULTIPLE  HOLOGRAPHIC  OPTICAL  ELEMENTS 

(submitted  to  Journal  of  the  Optical  Society  of  America  A,  1989) 

Abstract 

Using  the  stationary  phase  approximation,  holographic  optical 
elements  can  be  designed  which  perform  one-to-one  coordinate 
transformations.  All  coordinate  transformations  which  are  sufficiently 
differentiable  can  be  decomposed  into  two  transformations  in  series  and 
therefore  performed  using  two  holographic  elements.  Although  the 
method  presented  here  for  designing  the  two  required  holographic 
elements  is  not  explicit  for  the  general  case,  a  useful  necessary 
condition  for  the  decomposition  Is  given.  A  simple  explicit  solution 
based  on  separability  of  the  necessary  equations  for  the  first  of  the 
two  serial  transformations  Is  useful  In  some  cases. 

1.  Introduction 

Optical  processors  often  operate  on  data  that  are  presented  in  two 
spatial  dimensions.  In  some  applications,  the  data  need  to  be 
spatially  rearranged  at  some  point  In  the  processing.  These  data 
rearrangements  may  be  thought  of  as  one-to-one  coordinate 
transformations  (also  called  geometric  distortions  or  mappings). 
Conventional  optics  cannot  easily  perform  complicated  transformations, 
so  holographic  optical  elements  (HOEs)  are  often  used  for  coordinate 
transformations.  In  particular,  computer-originated  HOEs  are  used 
because  of  the  greater  ease  with  which  complicated  phase  functions  may 
be  fabricated.1,2 
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One  method  for  designing  the  HOEs  required  for  a  given  coordinate 
transformation  uses  the  stationary  phase  approximation  and  was 
developed  by  Bryngdahl.  To  review  this  method,  consider  the  optical 
system  shown  in  Fig.  1  where  the  Input  field  at  P1  in  the  front  focal 
plane  of  a  lens  L  Is  unlphase  and  has  an  amplitude  distribution  a^x, 
y).  A  HOE  with  phase  0(x,  y)  is  placed  in  plane  Pj.  The  complex 
amplitude  distribution  a2(s,  t)  found  in  the  back  focal  plane  P2  of  the 
lens  is  given  approximately  by  the  Fourier  transform: 

00  00 

a2(s'  t)  *  -  J  j  ajU,  y)  exp[lf(x,  y)  -  (xs  +  yt)]  dxdy  (1) 
-00  -00 

where  X  is  the  wavelength  and  f  is  the  lens  focal  length.  Following 
Bryngdahl,  this  integral  can  be  approximated  by  the  method  of 
stationary  phase5  to  show  that  a  point  (x,  y)  in  the  input  is  mapped  to 
a  point  (s,  t)  in  the  output  where 


Xf 

9#(x.i-y)l 

X 

s' 

2* 

9x 

■y. 

.t. 

Xf 

,2r 

9y 

The  transformation  Is  one-to-one  If  no  two  points  (Xj,  yj)  and  (x2,  y2) 
are  mapped  to  the  same  point  (s,  t) .  We  will  consider  only  one-to-one 
transformations. 

If  a  desired  coordinate  transformation  is 


x' 

s' 

g(x,  y)' 

y- 

*▼ 

.t. 

h(x,  y). 
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then  the  required  HOE  phase  ^(x,  y)  Is  found  by  solving  the  system  of 
partial  differential  equations: 


Xf 

2r  9x 


=  g(x,  y) 


hi  £4  =  h(x  v) 
2 r  ay  1  '  yj 


(4) 


It  will  be  assumed  throughout  that  the  region  in  the  x-y  plane  to  be 
transformed  Is  simply  connected  (1 . e . ,  that  It  has  no  holes)  and  that 
the  functions  g  and  h  are  continuously  differentiable  (i.e.,  that  the 
first  partial  derivatives  of  g  and  h  are  continuous).  Many  coordinate 
transformations  of  interest  satisfy  these  assumptions.1,2,3,4  Under 
these  assumptions,  the  system  of  Eq.  (4)  has  a  solution  if  and  only  if6 


a*(x.  y)  =  i-ll 

9ydx  5x3y 


which  is  equivalent  to 


9g.(x,  y)  =  3h(x,  y) 

9y  ax  ’  w 

A  coordinate  transformation  is  said  to  be  realizable  with  a  single  HOE 
designed  using  the  stationary  phase  approximation  if  and  only  if  Eq.  (6) 
is  satisfied. 

Unfortunately,  many  coordinate  transformations  of  interest  in 
1  2 

optical  processing  '  are  not  realizable  with  a  single  HOE.  An  example 

is  polar  formatting  for  image  formation  from  spotlight  mode  SAR  data. 

However,  two  HOEs  in  series  can  be  used  to  perform  this  coordinate 
2 

transformation.  In  this  case,  the  desired  transformation  is  decomposed 
into  two  serial  transformations,  each  of  which  is  realizable  with  a 
single  HOE.  [Note  that,  as  described  in  Ref.  2,  phase  correction  is 
performed  after  the  first  transformation  so  that  the  input  to  the  next 
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transformation  Is  unlphase  for  a  unlphase  Input.  The  phase  correction 
of  the  output  of  the  first  transformation  can  be  Included  In  the  HOE 
which  performs  the  second  transformation.  If  the  output  of  the  second 
transformation  must  be  unlphase  (for  example,  for  further  coherent 
optical  processing),  then  a  third  HOE  must  be  used  to  correct  the  output 
phase.] 

Given  this  example,  two  questions  naturally  arise:  First,  can  all 
one-to-one  coordinate  transformations  satisfying  the  above  assumptions 
of  a  simply  connected  region  and  continuous  differentiability  be 
achieved  by  the  use  of  multiple  HOEs  In  series?  Second,  If  so,  what  is 
the  minimum  number  of  HOEs  required?  In  this  paper,  these  two  questions 
are  investigated  and  answered.  In  Section  2,  the  questions  are 
developed  mathematically  and  expressed  as  a  system  of  partial 
differential  equations.  In  Section  3,  it  Is  shown  that  all  such 
transformations  can  be  achieved  with  at  most  two  HOEs  In  series. 
Section  4  discusses  orientation  preserving  and  reversing  transformations 
and  their  relation  to  the  earlier  results.  Section  5  applies  these 
results  to  an  example.  Conclusions  and  suggestions  for  further  research 
are  presented  In  Section  6. 

2.  Mathematical  Development 

We  begin  by  Investigating  the  conditions  under  which  a 
transformation: 


9s(x,  y)  .  3t(x,  y) 
3y  *  3x 


(8) 
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i.e.,  which  Is  not  realizable  with  a  single  HOE,  can  be  performed  by  two 
transformations, 


and 


x‘ 

u' 

u(x,  y)' 

■y. 

.V. 

v(x,  y). 

(9) 


u' 

s' 

p(u,  v)' 

-► 

11 

.V. 

,t. 

q(u,  v). 

(10) 


which  satisfy 


Mx,  *1  =  9v(x,_y) 
ay  8x 


(11) 


3p(u.  v)  =  9q(P^-V) 

dv  du 


(12) 


I.e.,  each  of  which  Is  realizable,  and  for  which,  of  course, 

s(x,  y)  =  p[u(x,  y),  v(x,  y)] 
t(x,  y)  =  q[u(x,  y),  v(x,  y)]  . 


(13) 


If  this  can  be  done,  then  the  transformation  can  be  performed  optically 
by  two  HOEs  In  series. 


We  begin  by  observing  that  a  one-to-one,  continuously  differentiable 
first  transformation  [given  by  Eq.  (9)]  onto  a  simply  connected  region 
can  always  be  chosen  that  is  realizable  [satisfies  Eq.  (11)]  and  has  a 
differentiable  inverse.  (The  differentiability  of  the  inverse  is  used 
later.)  Since  the  transformations  given  by  Eqs.  (7)  and  (9)  are  one-to- 
one,  then  a  second  one-to-one  transformation  [given  by  Eq.  (10)]  which 
satisfies  Eq.  (13)  will  always  be  well-defined.  The  question  then 
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becomes:  under  what  conditions  will  this  second  transformation  be 
realizable  [satisfy  Eq.  (12)]?  Clearly,  for  some  choices  of  the  first 
transformation,  the  required  second  transformation  will  not  be 
realizable.  However,  we  now  derive  a  system  of  partial  differential 
equations  which,  If  satisfied  by  the  first  transformation,  guarantees 
that  the  second  transformation  will  be  realizable. 

By  applying  the  chain  rule  for  partial  derivatives  to  Eq.  (13),  we 
arrive  at  a  system  of  equations  which  can  be  written  as: 


9s  9s' 
3x  3y 

at  at 
ax  ay. 


9p 

3u  dv 

93  M 
,3u  3v. 


3u  3u~ 
3x  9y 

3v  3v 
.3x  3y. 


(14) 


The  individual  matrices  have  the  same  elements  as  the  Jacobian 
determinants  for  each  transformation7  and  the  equation  expresses  the 
fact  that  the  Jacobian  matrix  of  the  total  transformation  is  the  product 
of  the  Jacobian  matrices  of  the  individual  transformations.  It  is 
Interesting  to  note  that  the  realizability  conditions  given  by  Eqs.  (11) 
and  (12)  are  equivalent  to  requiring  that  the  matrices  arising  from  the 
transformations  defined  In  Eqs.  (9)  and  (10)  be  symmetric. 

Our  next  step  Is  to  solve  Eq.  (14)  for  3p/3v  and  3q/3u  for 
substitution  Into  Eq.  (12).  This  can  be  done  by  inverting  the  matrix 
Eq.  (14),  i.e.,  multiplying  Eq.  (14)  from  the  right  by  the  inverse  of 
the  right-hand  matrix,  performing  the  matrix  multiplications,  and 
choosing  the  desired  terms.  The  result  Is: 


D  9jd  _  9s  9u  9s  3u 

9v  "  9y  9x  9x  3y 

D  =  9t  _  at  3v 
v  3u  ”  9x  9y  3y  9x 


(15) 


80 


ERIM 


where 


n  9u  9v  9u  9v 
u  '  3x  3y  "  3y  3x 


(16) 


D  is  the  Jacobian  determinant  for  the  first  transformation.  We  now  use 
the  assumed  differentiability  of  the  inverse  of  the  first 
transformation.  Since  D  is  never  equal  to  zero  for  a  transformation 

O 

with  a  differentiable  inverse,  Eq.  (15)  and  similar  expressions  for 
3p/9u  and  3q/9v  show  that  the  second  transformation  is  continuously 
differentiable. 


We  can  now  state  our  first  important  result:  For  the  total 

transformation  to  be  performed  as  a  series  of  two  realizable 
transformations  [i.e.,  for  Eqs.  (11),  (12),  and  (13)  to  be  satisfied],  a 
necessary  condition  is  that  the  first  transformation  in  the  series 
satisfies  the  system  of  partial  differential  equations: 

9u  .  6v 
9y  3x 

(17) 

at  av  at  av  t  as  au  as  au 

ax  ay  ay  ax  ay  ax  ax  ay  * 

The  first  equation  is  just  Eq.  (11)  and  the  second  is  obtained  from  Eq. 
(12)  by  using  Eq.  (15)  and  the  fact  that  D  is  never  equal  to  zero.  In 
addition,  the  first  transformation  must  be  one-to-one,  continuously 
differentiable,  onto  a  simply  connected  region,  and  have  a 
differentiable  inverse.  Together,  the  two  transformations  iri  series 
make  it  possible  to  perform  with  two  HOEs,  under  the  stationary  phase 
approximation,  the  total  transformation  that  was  not  realizable  with  a 
single  HOE. 
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3.  Existence  of  Solutions  to  the  Partial  Differential  Equation  System 

The  usefulness  of  the  result  of  Section  2  depends  on  the  existence, 
for  general  transformations  as  specified  by  Eq.  (7),  of  invertible 
solutions  to  the  system  of  partial  differential  equations  (PDEs)  given 
in  Eq.  (17).  Previous  work  has  shown  that  solutions  always  exist  to 

g 

linear  hyperbolic  PDE  systems  in  two  independent  variables.  We 
therefore  find  a  condition  for  which  the  PDE  system  is  hyperbolic.  We 
do  not  explicitly  find  a  general  solution  to  the  system  of  Eq.  (17). 

In  canonical  form,  Eq.  (17)  becomes:10 

9u  9v  n 
3y  ‘  9x  "  u 

(18) 

av  rati-1/  as  au  ras  ati  m 

ay  LaxJ  r  ay  ax  Lax  "  ayj  ax/  ~  u  * 

The  classification  of  the  PDE  system  is  determined  by  the  eigenvalues 
X(x,  y)  which  are  the  roots  of  the  characteristic  equation: 


Solving  Eq.  (19)  gives: 


(19) 
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The  PDE  system  will  be  hyperbolic  If  the  eigenvalues  are  real  and 
distinct  for  all  x  and  y  In  the  region  of  the  coordinate  transformation. 
This  is  the  case  when  the  expression  within  the  square  root  in  Eq.  (20) 
Is  greater  than  zero.  This  expression  can  be  manipulated  to  give: 

ras  at' ) 2  .  as  at 

lax  '  ayj  4  ay  ax  " 

(21) 

CM  +  Ml 2  +  a  f 8s  M  as  at] 
lax  ayj  4  lay  ax  "  ax  ayj  • 


We  now  need  to  require  only  that  the  second  term  be  greater  than  zero. 


A  sufficient  condition  for 

real  and 

distinct  eigenvalues  is  therefore: 

as 

3s 

3x 

3y 

<  0  . 

(22) 

at 

at 

ax 

ay 

The  condition  of  Eq. 

(22)  is 

very  important  because, 

if  it  is 

satisfied,  then  the  PDE  system  is 

hyperbolic.  For  a  hyperbolic  system 

in  two  independent  variables, 

it 

is  known  that  an  infinite 

family  of 

solutions  always  exists. 

This 

is 

proven  by  developing 

a  series 

solution  and  showing  that  It  converges  uniformly.  We  assume  that,  for 
the  one-to-one,  continuously  differentiable  coordinate  transformations 
onto  simply  connected  regions  which  are  of  Interest  here,  at  least  one 
member  of  this  family  will  have  an  inverse  that  is  also  differentiable. 

In  summary,  we  have  shown  that  the  total  transformation  can  be 
decomposed  into  two  transformations  in  series  if  a  solution  to  a  related 
PDE  system  [Eq.  (18)]  exists.  This  PDE  system  will  be  hyperbolic,  and 
therefore  an  infinite  family  of  solutions  will  exist,  if  the  total 
transformation  [defined  in  Eq.  (7)]  has  a  Jacobian  determinant  which  is 
negative  everywhere  in  the  region  to  be  transformed.  (It  will  be  shown 
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In  the  next  section  that  this  result  can  be  trivially  extended  to  total 
transformations  for  which  the  determinant  is  always  positive.) 

It  should  be  noted  that,  If  9t/9x  =  0,  then  this  approach  cannot  be 
used.  (This  means  that  t  Is  a  function  of  y  only  and  that  data  on  any 
line  parallel  to  the  x-axis  will,  after  transformation,  lie  on  a  line 
parallel  to  the  s-axis.)  However,  if  9s/9y  t  0,  then  a  second  approach 
is  possible  which  begins  by  dividing  Eq.  (17)  by  9s/9y.  If  both  9t/9x 
and  9s/9y  =  0,  then  the  desired  transformation  satisfies  the 
realizability  condition  and  there  is  no  need  to  decompose  it  into  two 
serial  transformations. 

4.  Orientation  Preserving  and  Reversing  Transformations 

To  understand  the  results  of  Section  3,  we  need  to  define 
orientation  preserving  and  orientation  reversing  transformations.  A 
transformation  is  said  to  be  orientation  preserving  if  a  point 
traversing  a  closed  curve  clockwise  in  the  original  coordinate  system 
also  traverses  the  transformed  curve  clockwise.11  A  transformation  is 
orientation  reversing  if  the  point  traverses  the  transformed  curve 
counterclockwise.  Orientation  preserving  transformations  are  like 
ordinary  images  while  orientation  reversing  transformations  are  like 
mirror  images. 

It  is  also  true  that  a  transformation  is  orientation  preserving  if 
the  Jacobian  determinant  is  positive  and  orientation  reversing  when  the 
determinant  is  negative.11  By  Eq.  (22),  we  require  the  determinant  to 
be  negative  for  the  resulting  POE  system  to  be  hyperbolic  and  therefore 
have  a  solution,  so  we  are  in  effect  requiring  that  the  transformation 
be  orientation  reversing. 

This  seeming  restriction  to  transformations  which  are  orientation 
reversing  can  easily  be  relaxed.  An  oriention  preserving  transformation 
can  be  made  to  be  orientation  reversing  simply  by  reversing  the 
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direction  of  one  of  the  output  coordinates  (thus  making  it  a  mirror 
image).  In  an  optical  processor,  this  usually  means  just  redefining  the 
output  data  coordinates.  If,  for  some  reason,  the  coordinates  cannot  be 
redefined,  then  a  third  coordinate  transformation  which  turns  the 
orientation  preserving  transformation  into  an  orientation  reversing  one 
is  necessary.  One  possible  transform  is: 

(23) 

which,  of  course,  is  realizable  with  a  single  HOE.  It  should  also  be 
noted  that  the  requirement  for  an  orientation  reversing  transformation 
has  only  been  shown  to  be  a  sufficient,  not  a  necessary,  condition  for 
the  existence  of  a  solution  to  the  PDE  system. 

5.  Illustrative  Example 

To  better  understand  the  above  results,  it  may  be  helpful  to 

consider  a  specific  example.  For  polar  formatting,  the  desired 

2 

coordinate  transformation  is: 


x' 

s' 

x  cos  y' 

= 

y. 

.t. 

,x  sin  y. 

(24) 


where  x-y  and  s-t  are  both  orthogonal  coordinates.  This  transformation 
performs  polar  formatting  because,  if  we  identify  the  r  and  6  of  polar 
coordinates  with  x  and  y,  respectively,  then  r  and  9  are  rectilinear 
coordinates  in  the  x-y  system,  but  polar  coordinates  in  the  s-t  system. 
The  transformation  functions  are: 

s(x,  y)  =  x  cos  y 

(25) 

t(x,  y)  =  x  sin  y  . 
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The  region  of  transformation  has  to  be  restricted  to  x  >  0  or  x  <  0  for 
the  transformation  to  be  one-to-one. 

Upon  presentation  of  a  coordinate  transformation  such  as  this,  the 
first  step  Is  to  determine  whether  or  not  Eq.  (6)  is  satisfied.  The 
partial  derivatives, 


-x  sin  y 


3t(x.  y) 
3x 


sin  y  , 


are  not  equal  and,  hence,  the  transformation  is  not  realizable  with  a 
single  HOE. 

The  second  step  Is  to  determine  if  It  is  orientation  reversing 
(i.e.,  If  Eq.  (22)  is  satisfied).  Again,  it  is  easily  shown  that  it  is 
not.  Therefore,  to  make  use  of  the  results  of  Sections  2  and  3,  the 
desired  transformation  functions  can  be  redefined  as: 

s (x ,  y)  =  x  cos  y 

(27) 

t(x,  y)  =  -x  sin  y  . 

(Other  choices  are  also  possible.)  The  minus  sign  in  Eq.  (27)  is  used 
so  that  Eq.  (22)  Is  satisfied,  and  the  transformation  is  orientation 
reversing,  for  x  >  0.  The  only  physical  effect  of  the  minus  sign  is  to 
invert  one  of  the  axes  In  the  output  plane. 

Using  Eq.  (25),  the  canonical  form  [Eq.  (18)]  for  this  example  is: 


3u  3v  n 
3y  "  3x  '  u 


Sy  +  I’  Sin  yl  1(X  s1n  y  S  +  [(X  +  !)  cos  yl  fx)  = 


(28a) 


(28b) 
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A  solution  of  Eq.  (28)  can  be  found  by  assuming  that  u  =  g(x)  and  v  = 
h(y).  Then  Eq.  (28a)  is  satisfied  trivially  and  Eq.  (28b)  becomes: 


dh 

dy 


(29) 


Eq.  (29)  is  easily  solved  and  the  solutions  give  a  first 
transformation  of  the  form: 


x' 

u' 

c  In  x  +  d' 

-► 

= 

•yJ 

.V. 

cy  +  e 

(30) 


where  c,  d,  and  e  are  constants.  This  transformation  is  one-to-one, 
continuously  differentiable,  onto  a  simply  connected  region,  and  has  a 
differentiable  inverse.  Using  Eqs.  (27)  and  (30),  the  second 
transformation  is: 


u 

s' 

exp  u  cos  v' 

= 

•  Vj 

.t. 

-exp  u  sin  v. 

where  the  constants  have  been  determined  to  have  the  values  c  =  1,  d  =  e 

=  0.  Note  the  minus  sign  that  comes  from  Eq.  (27)  and  makes  the  total 

transformation  orientation  reversing.  The  polar  formatting 

transformation  has  therefore  been  separated  into  two  transformations 

performed  in  series  that  are  realizable  [i.e.,  that  satisfy  Eqs.  (11), 

(12),  and  (13)].  These  transformations  have  been  implemented  with  HOEs 

2 

and  performed  as  predicted. 

Note  that  this  example  was  solved  by  trying  a  first  coordinate 
transformation  which  was  separable,  i.e.,  u  was  a  function  of  x  only  and 
v  a  function  of  y  only.  This  method  may  be  useful  in  designing  HOEs  for 
other  transformations.  As  another  aid  to  design,  it  is  important  to 
note  that  a  system,  designed  to  perform  a  specific  coordinate 
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transformation,  will  perform  the  Inverse  transformation  if  the 
propagation  direction  of  the  light  is  reversed.  In  this  case,  the  phase 
correction  HOEs  discussed  in  Section  1  and  the  designed  coordinate 
transformation  HOEs  simply  reverse  their  roles. 

6.  Conclusion 

Because  of  the  need  to  perform  data  transformations  optically,  the 
general  problem  of  one-to-one  coordinate  transformations  realizable  with 
multiple  HOEs  was  investigated.  Subject  to  some  requirements  on  the 
transformation  that  can  be  met  in  most,  perhaps  all,  practical  cases, 
any  one-to-one,  continuously  differentiable  transformation  over  a  simply 
connected  region  can  be  performed  with  at  most  two  HOEs  designed  using 
the  stationary  phase  approximation.  One  design  method  was  described 
that  was  useful  for  the  polar  coordinate  transformation  and  may  be 
useful  in  other  cases  as  well.  However,  further  research  is  necessary 
to  develop  a  method  for  designing  the  two  HOEs  that  would  be  applicable 
in  all  cases.  Another  problem  of  interest  is  the  general  theory  of 
many-to-one  and  one-to-many  optical  transformations.  Of  possible 
mathematical  interest  would  be  the  development  of  a  proof  that  at  least 
one  member  of  the  infinite  family  of  solutions  to  a  hyperbolic  system  of 
PDEs  has  a  di fferentiable  inverse. 

This  research  was  supported  by  the  U.S.  Army  Research  Office. 
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Fig.  1.  Optical  system  for  coordinate  transformation.  The  lens  L  has 
focal  length  f  and  Fourier  transforms  plane  to  plane 
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ABSTRACT 


Computer-generated  holograms  (CGHs)  are  used  in  a  number  of  important  optical  technology  applica¬ 
tion  areas  such  as  holographic  optical  elements,  optical  processing  and  computing,  optical  testing,  image 
and  information  display,  beam  forming,  and  beam  scanning.  Many  different  CGH  fabrication  devices  (e.g., 
laser  beam  scanners,  electron  beam  writers)  and  facilities  have  been  developed  and  are  in  use.  However, 
none  of  these  devices  ideally  suit  the  requirements  (e.g.,  resolution,  space-bandwidth  product,  recording 
material)  of  many  CGH  applications.  Furthermore,  the  access  of  many  researchers  to  these  facilities  is 
limited  and  the  technical  support  available  is  often  poor.  New  facilities  specifically  designed  for  CGH  fab¬ 
rication  would  better  serve  the  needs  of  the  CGH  research  and  development  community.  The  requirements 
for  a  successful  CGH  fabrication  facility  including  appropriate  technical  support  for  users  are  established. 

1.  INTRODUCTION 


The  possibility  of  a  national  facility  for  the  fabrication  of  computer-generated  holograms  (CGHs)  has 
been  considered  for  the  past  decade.  Some  of  the  possible  benefits  of  a  national  CGH  facility  are:  (1)  it 
would  fabricate  high  quality  CGHs  for  use  by  the  research,  development,  and  manufacturing  communities, 
(2)  it  would  avoid  further  expensive,  needless  duplication  of  facilities,  (3)  it  would  speed  the  development 
and  application  of  diffractive  optics  techniques  and  devices,  and  (4)  it  would  improve  system  performance 
in  several  important  application  areas  (e.g.,  holographic  optical  elements,  optical  interconnect  devices  for 
high  speed  parallel  processors,  and  invariant  correlation  Filters  for  object  detection  and  recognition). 

Although  many  facilities  are  currently  fabricating  CGHs  of  various  types,  all  of  these  facilities  suffer 
from  one  or  more  of  the  following  deficiencies:  (1)  performance  limitations  (e.g.,  resolution,  space- 
bandwidth  product,  or  recording  speed),  (2)  cost,  (3)  CGH  encoding  or  materials  restrictions  (e.g.,  restric¬ 
tion  to  binary  amplitude  or  phase,  or  restriction  to  recording  fringes),  (4)  data  format  inefficiencies,  and 
(5)  lack  of  technical  support  in  CGH  design  and  optics  in  general.  Clearly,  no  existing  facility  can  yet  be 
considered  to  be  acceptable  as  a  national  CGH  fabrication  facility. 

This  paper  makes  an  effort  to  take  the  next  step  toward  a  better,  and  possibly  national,  CGH  fabrica¬ 
tion  facility.  A  broad  overview  of  CGH  applications  and  technology  is  given  in  Section  2.  CGH  fabrica¬ 
tion  methods  and  devices  are  discussed  in  Section  3.  Finally,  in  Section  4,  the  requirements  which  a  suc¬ 
cessful  CGH  fabrication  facility  must  meet  are  given  and  an  in-depth  study  of  user  requirements  is 
recommended. 

2.  REVIEW  OF  CGH  TECHNOLOGY 


The  range  of  applications  which  can  potentially  take  advantage  of  CGH  is  very  wide.  Some  of  the 
applications  are  shown  in  Figure  1.  They  have  been  grouped  according  to  whether  the  CGH  is  used  to 
modify  only  the  phase  of  a  wavefront  or  to  modify  both  amplitude  and  phase. 

CGH  can  be  classified  in  many  ways.  Two  possible  methods  are  by  material  type  and  by  type  of  dif¬ 
fraction  effect  utilized.  The  different  types  of  CGH  under  each  of  these  categories  are  shown  in  Figure  2. 
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CGH  APPLICATIONS 
Phase-only  wavefront  modification 


Conventional  optics 

HOE  •  lens,  mirrors,  phase  correctors 
Testing  of  aspherics 

Wavefront  creation  for  COHOE  fabrication 
Scanners 

Heada-ap/helmet  displays 
Light  redistribution 

Beam  intensity  profile  shaping 
Data  formatting  -  coordinate  transformation 
Digital  data  processor  interconnects 
Beam  combiners,  dividers 

Amplitude  and  phase  wavefront  modification 


Displays 

2-D,  3-D  images 
Sighting  devices,  reticles 
Spatial  filters  for  pattern  recognition 
Data  memory,  storage 


Figure  1.  The  wide  range  of  CGH  applications. 


CGH  CLASSIFICATION 

Material  types 

Amplitude-only 

Continuous 

Binary 

Phase-only 

Continuous 

Binary 

Amplitude  and  phase 

Diffraction  effects 
Scalar 

Planar  media  and  surface  relief 
Volume  Bragg  diffraction 
Vector 


Figure  2.  Examples  of  CGH  cli  uification  methods. 
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CGH  can  b«  copied  in  several  ways.  A  contact  copy  can  be  made,  usually  for  the  purpose  of  trans- 
f erring  the  CGH  onto  a  new,  more  suitable  material.  Another  possibility  is  the  COHOE  or  Computer- 
Originated  Holographic  Optical  Element.  In  this  copying  method,  the  wavefront  reconstructed  by  the  CGH 
is  interfered  with  a  second  wavefront.  For  example,  if  this  second  wavefront  is  planar,  then  the  effect  is 
simply  to  copy  the  phase  of  the  CGH  while  adding  a  high  carrier  frequency  to  achieve  higher  diffraction 
efficiency  or  a  larger  diffraction  angle.  Another  examp' e  is  for  the  second  wavefront  to  be  spherical  and 
of  high  numerical  aperture  and  the  CGH  wavefront  to  represent  a  desired  phase  correction;  then  the 
resulting  COHOE  can  be  a  high  numerical  aperture  aspheric.  In  the  latter  example,  the  space-bandwidth 
product  required  of  the  CGH  is  reduced  by  using  a  nonplanar  wavefront  generated  by  conventional  optics. 

Because  of  the  wide  range  of  applications  in  which  CGHs  may  be  used,  there  is  also  a  wide  range  in 
the  performance  required  of  CGHs  (see  Figure  3).  The  specifications  shown  in  the  low  range  might  apply 
to  a  CGH  for  a  Fourier- transform-domain  matched  filter  in  a  coherent  optical  correlator  and  those  in  the 
high  range  to  an  aspheric  holographic  combiner  for  a  head-up  display.  For  many  applications,  the 
required  signal-to-noise  ratio  has  not  been  adequately  analysed,  hence  the  question  marks  in  the  table  rath¬ 
er  than  numerical  values.  In  addition,  the  diffraction  angle  required  of  a  CGH  can  vary  from  sero  to 
nearly  180*  depending  on  the  application.  Although  CGHs  have  been  most  commonly  used  in  the  wave¬ 
length  region  from  the  ultraviolet  to  the  infrared,  CGHs  have  also  been  made  and  used  at  x-ray  and  mi¬ 
crowave  wavelengths.  It  is  neither  expected  or  required  that  any  one  CGH  fabrication  device  would  meet 
all  the  specifications  shown  in  Figure  3. 


RANGE  OF 

CGH  PERFORMANCE  SPECIFICATIONS 

Ranee 

Low  High 

Resolution 

100  fim 

0.05  /im 

Space-bandwidth  product 

102  x  102 

10®  x  108 

Sise 

(1  mm)2 

(50  cm)2 

Diffraction  effrciency 

<1% 

~100% 

Phase  accuracy 

X 

X/50 

Signal-to-noise  ratio 

7 

7 

Figure  3.  The  wide 

range  of  CGH  specifications 

for  different  applications. 

There  are  also  a  wide  variety  of  methods  for  encoding  the  desired  two-dimensional  distribution  of  am¬ 
plitude  and  phase  in  a  CGH.  The  major  encoding  methods  are  shown  in  Figure  4.  Each  me -hod  has  its 
own  advantages  and  disadvantages.  For  the  case  of  CGH  used  for  two-dimensional  image  reconstruction, 
the  inaccuracies  inherent  in  most  of  these  methods  have  been  analysed. 
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•  MAJOR  CGH  ENCODING  METHODS 

Amplitude  and  phase  directly  •  ROACH 
»  Phase  directly  -  Kinoform 

Continuous 

Binary 

Phase  on  carrier 

Amplitude  transmittance 
Fringes 

Continuous  -  Burch 
Binary  •  Lee 
Detour-phase 

Continuous  •  Burckhardt  and  Lee 
Binary  •  Lohmann 
Phase  transmittance 

Phase  versions  of  above 


Figure  4.  A  classification  of  the  major  CGH  encoding  methods. 


As  has  been  the  case  in  each  CGH  characteristic  examined  in  this  section,  the  range  of  CGH  recording 
media  is  also  very  wide.  Some  of  the  principal  types  are  listed  in  Figure  5.  It  is  important  to  note 
that,  even  when  the  CGH  can  not  be  recorded  in  a  material  suitable  for  use  in  the  application,  it  can 
usually  be  replicated  in  a  suitable  material.  Not  included  in  this  list  are  the  various  real-time  spatial 
light  modulators  that  can  be  used  to  record  CGHs.  They  are  ordinarily  restricted  to  uses  in  which  the 
hologram  must  be  recorded  in  near  real-time,  and  are  therefore  not  appropriate  for  a  CGH  facility.  (They 
also  typically  have  modest  performance  specifications). 


CGH  RECORDING  MEDIA 

Photographic  emulsions 
Amplitude 
Bleached 
Gelatin 

Silver  halide 
Dichromated 

Photoresist,  electron  resist 
Photopolymers 
Multiemulsion  -  ROACH 
Machinable  materials 

Replication  materials 
Above  plus: 

Glass,  quarts 
Si,  GaAs,  Ge 
Metal 


Figure  S.  Examples  of  CGH  recording  and  replication  media. 
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3.  CGH  FABRICATION  DEVICES 


A  large  number  of  devices  have  been  used  to  fabricate  CGH.  The  principle  methods  and  devices  are 
listed  in  Figure  0.  In  addition  to  the  categories  listed  in  Figure  0,  the  devices  can  be  divided  into  two 
classes:  (1)  those  which  were  built  for  another  purpose  and  were  adapted  for  CGH  fabrication  and 
(2)  those  which  have  been  designed  and  built  specifically  for  CGH  fabrication.  Plotters,  printers,  rotating 
drum  scanners,  some  translating  flat  bed  scanners,  optical  pattern  generators  for  microelectronic  mask 
making,  electron  beam  writers,  and  computer-controlled  milling  machines  are  members  of  the  first  class. 
Although  useful  CGH  have  been  made  with  all  these  devices,  better  performance  (e.g.,  smaller  writing  spot 
sise,  greater  space-bandwidth  product,  more  appropriate  data  formatting  and  CGH  encoding)  could  be 
obtained  if  the  device  were  redesigned  specifically  for  CGH  fabrication.  Galvanometer  scanners,  devices 
which  image  a  CRT  to  the  recording  medium,  and  devices  which  interfere  two  computer-controlled 
wavefronts  belong  to  the  second  class.  Although  some  of  these  devices  are  excellent  for  fabricating  specific 
types  of  holograms,  they  all  have  performance  limitations  which  make  them  unsuitable  for  fabricating  a 
wide  variety  of  CGH  types. 


CGH  FABRICATION  METHODS  AND  DEVICES 

Drawing 

Plotter 

Printer 

Optical  writing 

With  one  beam 

Galvanometer 
Rotating  drum 

Translating  flat  bed  •  microdensitometer 
Optical  pattern  generator 
With  two  beams 

Interference  of  two  wavefronts 
With  CRT  image 
Electron  beam  writing 
Machining  of  surface  profile 

Figure  0.  The  major  CGH  fabrication  methods  and  representative  device  types. 


For  example,  consider  devices  which  use  a  laser  beam  to  write  on  a  photosensitive  film  placed  on  a 
rotating  drum.  These  devices  have  the  advantages  of  (1)  continuous  modulation  of  the  intensity  of  the  la¬ 
ser  beam  allowing  a  continuous  amplitude  recording,  (2)  being  relatively  inexpensive  to  purchase  and 
maintain,  (3)  fast  recording  speed,  (4)  potentially  high  space-bandwidth  product,  and  (5)  moderate  resolu¬ 
tion.  Disadvantages  are  (1)  positional  inaccuracy  can  cause  phase  errors,  (2)  current  use  with 
photographic  film  produces  a  CGH  with  (a)  phase  errors  due  to  thickness  variations  (correctable  by 
mounting  in  an  index  matching  fluid)  and  (b)  low  diffraction  efficiency  (correctable  by  copying  onto  a 
high  efficiency  phase  material),  and  (3)  restriction  to  a  raster  exposure  pattern  which  is  not  optimum  for 
writing  CGH  which  consist  of  curved  fringe  patterns. 

The  currently  popular  electron  beam  recorder  which  is  used  for  mask  generation  for  microelectronics  is 
another  example.  It  has  the  advantages  of  (1)  high  (submicron)  resolution,  (2)  high  space-bandwidth 
product,  and  (3)  allowing  optics  to  build  on  microelectronics  technology.  However,  current  use  of  this  de¬ 
vice  suffers  from  several  disadvantages:  (1)  the  recording  is  principally  binary,  (2)  the  devices  and  their 
maintenance  are  expensive,  (3)  recording  speed  is  only  moderate  to  low,  and  (4)  much  software  develop¬ 
ment  is  need  to  allow  easier  user  access,  to  permit  efficient  data  formating,  and  to  make  use  of  the 
random-access  ability  to  write  curved  CGH  fringe  patterns. 
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4.  CQH  FACILITY  REQUIREMENTS 

The  brief  review  of  CGH  applications  and  technology  given  in  Section  2  indicates  that  multiple  issues 
are  involved  in  the  choice  of  a  CGH  fabrication  method  and  in  the  design  of  a  CGH  fabrication  facility. 
As  indicated  by  the  examples  of  fabrication  devices  given  in  Section  3,  all  current  facilities  have  some 
difficiencies.  Designers  and  builders  of  future  facililites  should  consider  the  CGH  facility  requirements 
summarised  in  Figure  7.  First,  the  projected  facility  must  offer  high  performance  levels.  High  perfor- 
mance  is  not  sufficient  for  success  by  itself,  however.  The  facility  must  also  have  the  flexibility  to 
accomodate  the  various  encoding  methods  and  recording  materials  required  by  different  applications. 
Finally,  the  facility  hardware,  software,  and  technical  and  management  staff  must  be  user-oriented.  Only 
if  all  three  of  these  requirements  are  met  will  it  be  asked  to  produce  a  sufficient  volume  of  CGHs  to  jus¬ 
tify  the  financial  investment  in  a  new  CGH  fabrication  facility. 


REQUIREMENTS  FOR  A  SUCCESSFUL  CGH  FACILITY 

High  performance 

High  space-bandwidth  product 
High  resolution 
Low  errors 
Fast  recording 

Flexibility  to  accomodate  various 
Encoding  methods 
Recording  media 
User-oriented 

Flexible  CGH  description 
Low  cost 

Reasonable  delivery  time 
Technical  support 

Figure  7.  The  three-fold  requirements  for  a  successful  CGH  facility. 


It  is  not  possible  at  this  time  to  determine  which  fabrication  method  (or  methods),  hardware,  software, 
and  support  structure  should  be  chosen  for  a  new  CGH  fabrication  facility.  The  CGH  performance 
requirements  of  potential  users  are  complex  and  require  more  study.  When  this  has  been  accomplished,  a 
trade-off  study  can  be  performed  in  an  effort  to  select  a  CGH  fabrication  method  and  facility  that  best 
meets  user  community  needs. 
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